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This paper deals with multidimensional dynamic risk measures induced by conditional g- 
expectations. A notion of multidimensional (/-expectation is proposed to provide a multidimen- 
sional version of nonlinear expectations. By a technical result on explicit expressions for the 
comparison theorem, uniqueness theorem and viability on a rectangle of solutions to multidi- 
mensional backward stochastic differential equations, some necessary and sufficient conditions 
are given for the constancy, monotonicity, positivity, homogeneity and translatability properties 
of multidimensional conditional (/-expectations and multidimensional dynamic risk measures; we 
prove that a multidimensional dynamic (/-risk measure is nonincreasingly convex if and only if 
the generator g satisfies a quasi-monotone increasingly convex condition. A general dual repre- 
sentation is given for the multidimensional dynamic convex (/-risk measure in which the penalty 
term is expressed more precisely. It is shown that model uncertainty leads to the convexity of 
risk measures. As to applications, we show how this multidimensional approach can be applied 
to measure the insolvency risk of a firm with interacted subsidiaries; optimal risk sharing for 
7-tolerant (/-risk measures is investigated. Insurance (/-risk measure and other ways to induce 
(/-risk measures are also studied at the end of the paper. 

Keywords: multidimensional dynamic convex risk measure, backward stochastic differential 
equation, (/-expectation, insolvency risk, stochastic interaction, risk sharing. 

1. INTRODUCTION 

Quantifying the risk of the uncertainty in the future value of a portfolio is one of the key 
tasks of risk management. This quantification is usually achieved by modeling the uncertain 
payoff as a random variable, to which then a certain functional is applied. Such functionals are 
usually called risk measures. Coherent risk measures have been introduced by Artzner, Del- 
baen, Eber, and Heath (1997, 1999), and further developed by Delbaen (2000, 2002); sublinear 
risk measures by Frittelli (2000); convex risk measures by Follmer and Schied (2002a, b) and 
Frittelli and Rosazza Gianin (2002). In all these papers, static risk measures are under con- 
sideration. However practically we are interested in monitoring the riskiness of a position £ at 
any intermediate time t between the initial date and the final T. The dynamical risk mea- 
sure represents the riskiness at time t, taking into account all the information available up to 
date t. A class of dynamic risk measures, called g-risk measure, can be obtained by means of 
conditional (/-expectations (Jiang (2008); Rosazza Gianin (2006)). In particular, any dynamic 
g-risk measure satisfies the "time-consistency property". Note that the notion of time consis- 
tency first appeared in Peng (2005) under the Brownian filtration. See Kloppel and Schweizer 
(2007), Bion-Nadal (2008a, 2009a) and references therein for more recent papers dealing with 
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dynamic risk measures on a filtered probability space. Particularly Bion-Nadal (2008a, 2009a) 
characterized the time consistency by a "cocycle condition" for the minimal penalty function 
of a convex risk measure. And using the theory of BMO 1 martingales, she also provided a new 
class of time consistent risk measures even with jumps. A characterization of time consistency 
for a concave monetary utility function in terms of its acceptance set was presented in Cherid- 
ito, Delbaen, and Kupper (2006). Further discussions on dynamic risk measures can be found in 
Barrieu and El Karoui (2006), Frittelli and Rosazza Gianin (2004), Jobert and Rogers (2008), 
Roorda, Schumacher, and Engwerda (2005), among many others. The variety and complexity of 
the many risk factors and their interaction give rise to the need for a multidimensional approach 
to risk measurement. The nature of complex marketable assets has increased the requirements 
for necessary analytical methods. These methods must deal with the multiplicity of risk sources 
and their correlations. Many economists and mathematical finance scholars have attempted to 
achieve a multidimensional framework for risk analysis (Fong and Vasicek (1997); Jouini, Med- 
deb, and Touzi (2004); Kulikov (2007)). The notion of multidimensional static risk measure 
was introduced in Jouini, Meddeb, and Touzi (2004) and Kulikov (2007). The objective of the 
paper is to deal with multidimensional dynamic risk measures induced by multidimensional con- 
ditional g-expectations. To obtain a preliminary intuition of g-risk measures, we give a brief 
introduction. Let (Y , . . . ,Y n ) be the unique solution to multidimensional backward stochas- 
tic differential equation (BSDE for short) (2.1). (Y 1 ,...,Y n ) is also called multidimensional 
conditional ^-expectation here, denoted by (S^ H-Ft], • • • , £g n [-\J-'t]) ■ The (/-risk measure is then 
defined as (pf . . . , pf [£»]) = • • • , & n 9n [-fV*]) • If £, k ,k = 1, . . . n, represents 

the potential loss of the kth subsidiary of a firm, then pf k calculates its risk at time t. Due to 
the interaction among subsidiaries, pf 1 [£ ], . . . , pf n [£ n ] are related to each other through the gen- 
erator (gi, . . . g n ). See the application section for more intuitive examples. Note that the (7-risk 
measure depends heavily on its generator g. The coefficient g can be interpreted as infinitesimal 
risk measure over a time interval [t, t + dt] as pointed out in Barrieu and El Karoui (2006). The 
choice of g is crucial. Particularly, Chen and Epstein (2002) showed that the function g could 
depend on the preferences of investors. For one dimensional case, based on results in Coquet, 
Hu, Memin, and Peng (2002), Peng (2005), Delbaen, Peng, and Rosazza Gianin (2010), authors 
there showed that almost any dynamic risk measure under a Brownian filtration (satisfying cer- 
tain domination condition) can be represented by a conditional (/-expectation and an explicit 
representation of the penalty term of general dynamic concave utilities with continuity from 
above can be obtained by applying the theory of ^-expectation. This is also a motivation for us 
to investigate multidimensional risk measures induced by multidimensional ^-expectations. 

(/-expectations which were introduced in Peng (1997) provide an excellent explanation for the 
famous Allais's paradox, (/-expectation is an extension of the classical mathematical expectation. 
It preserves much of the properties of standard mathematical expectation except the linearity. 
The notion of (/-expectation can be applied to describe nonlinear recursive utilities, risk measures 
and the pricing of contingent claims, see El Karoui, Peng, and Quenez (1997) and references 
therein for an overview of applications. 

As far as now, few papers studied multidimensional (/-expectations. One of the important 
reasons is that we can not give an explicit necessary and sufficient condition for the comparison 
theorem of solutions of multidimensional backward stochastic differential equations, and the 
form of the multidimensional comparison theorem is more complicated than the 1-dimensional 



1 The space of functions of bounded mean oscillation. 
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one. Starting from the point of backward stochastic viability property (BSVP) in Buckdahn, 
Quincampoix, and Rascanu (2000), Hu and Peng (2006) gave a necessary and sufficient condition 
under which the comparison theorem holds for multidimensional BSDEs. This paper develops an 
explicit expression for the necessary and sufficient condition of the multidimensional comparison 
theorem which leads the multidimensional comparison theorem to be more applicable. 

Note that the multidimensional ^-expectation is not a notion just by putting together several 
one dimensional (^-expectations. It has its own complications due to the high dimensions. For 
a one dimensional g-risk measure induced by a ^-expectation, the monotonicity holds naturally 
due to the comparison theorem of BSDEs. But for multidimensional g-risk measures we need the 
"quasi-monotonicity condition" of the generator g throughout the paper. For the multidimen- 
sional case any axiom is a strong condition. In fact one should add axioms to multidimensional 
g-risk measures carefully. We will see later that if one expect to find examples of g- risk measures 
which satisfy the axioms of monotonicity and translation invariance, one may discover that the 
multidimensional g-risk measure is in fact n one dimensional g-risk measures put together. We 
extend notions of convex and coherent risk measures in this paper by dropping the axiom of 
translation invariance. 

We also present some applications of multidimensional g-risk measures. As interpreted in 
Kirman (1997) and Kirman, Oddou, and Weber (1986), economic agents do, in fact, communicate 
with each other and learn from each other. They also infer information from the actions of others 
and, most importantly, in most markets they trade with each other. Within a multi-business firm, 
subsidiaries interact with each other. If we understand how their mutual interactions influence 
the decisions, we can model their insolvency risk as a multidimensional BSDE. Recall that 
Jarrow (2002) measured the insolvency risk of a firm by the cost of a put option written on the 
firm's net value. Such an option guarantees that the value of the firm's assets will always be 
greater than or equal to the value of the firm's liabilities. So the firm could write a set of related 
put options for her subsidiaries to protect them against insolvency. The measure defined by the 
price of the put option is called "Put Premium Risk Measure". The acceptance set generated 
by this put premium risk measure generates a risk measure satisfying the axiom of translation 
invariance but may violating time consistency except that the original put premium risk measure 
also satisfies the axiom of translation invariance. This risk measure characterizes the minimal 
capital required to protect the firm's subsidiaries against insolvency. 

To deduce properties of multidimensional g-risk measures, we have to make some mathemat- 
ical preliminaries first. Readers who are interested only in financial applications can jump to 
section 4 and 5 with little knowledge about BSDEs. The paper is organized as follows. Some 
results on backward stochastic viability property are recalled in section 2. Section 3 contains 
the explicit expressions for the multidimensional comparison theorem and viability on a rect- 
angle. Some basic properties and Jensen's inequality for multidimensional (^-expectations are 
studied in section 4. Section 5 deals with multidimensional dynamic risk measures induced by 
multidimensional conditional g-expectations. The multidimensional g-risk measure is shown to 
be a multidimensional dynamic time consistent risk measure. We prove that a multidimen- 
sional dynamic (7-risk measure is nonincreasingly convex if and only if the generator g satisfies 
a quasi-monotone increasingly convex condition. A general dual representation is also given for 
a multidimensional dynamic convex (7-risk measure in which the penalty term is expressed more 
precisely. An application of our multidimensional approach is taken into account to measure the 
insolvency risk of a multi-business firm with interacted business units. Optimal risk sharing for 
7-tolerant g-risk measures is investigated. Insurance g-risk measure, multidimensional g-pricing 
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mechanism and stochastic differential utility are also considered. 



2. BACKWARD STOCHASTIC VIABILITY PROPERTY AND 
MULTIDIMENSIONAL COMPARISON 

Let (Bt)t£[o : T] be a standard d-dimensional Brownian motion on a probability space (£l,F, P) 
and (^ r t)tg[o,T] be the augmented Brownian filtration generated by (-Bt)ig[o,T]- Let T < be a 
fixed time. For any Euclidean space H, we define for p > 1, Lj-(0, T; i7):={all -ff-valued and 

J^-adapted stochastic processes such that E \(ft\ p dt < oo}; 5j-(0, T; H):={ail continuous 

processes in Ejr(0,T; H) such that E [sup 0<t<T \tpt\ p dt\ < oo}; LP {Ft; H):={the space of all 
.ff-valued and J^-measurable random variables such that E'[|£| p ] < oo}. We say x > y, for 
x, y 6 i? n , if Xj > yi, for all £ = 1,2,..., n. Let e^, i = l,...,n, denote the orthogonal basis of 

vectors in R n . Va 6 i? n , c G R nxd , we use norms |a| = V^iLi |°i| 2 ) l c l = y^Li X^j=i l c «jl 2 - 
Consider the following n-dimensional BSDE: 



(2.1) Vt = t + g{s,y s ,z s )ds- z s dB s ,0<t<T, 

Jt Jt 

where the terminal variable £ G L 2 {Ft; i? n ) and the function g : Q x [0,T] x R n x R nxd i — >. R n 
is progressively measurable. We call (£,g,T) the standard parameters of BSDE (2.1). 

For the function g, we assume throughout the paper that 

(HI) P-a.s., V{y,z), t — > g(t,y,z) is continuous; 

(H2) There exists a constant L ^ such that P-a.s., we have: Vt, V(y, y'), V(z,z / ), 

Iff (*, y,z) - g (t, y , z')\< L [\y - y'\ + \z- z'\) ; 
(H3) sup < t < T 0,0)| G L?(F T ;R n ). 

Let K be a nonempty closed convex set of R n . We define Hjc{a) as the projection of point a 
onto K, and dx(-) the distance function of X. The following proposition comes from Buckdahn, 
Quincampoix, and Rascanu (2000). Throughout this paper, it is understood that an equation 
or an inequality holds true always in the sense of P-a.s.. 



Proposition 2.1 Suppose that g satisfies (HI) ~(H3). Then the following are equivalent: 

(i) For any t G [0,T], V£ G L 2 {F t ;K), the unique solution (Y, Z) G 5jr(0, t; R n ) xL 2 r (0, t; R nxd ) 
to the following BSDE over time interval [0,t]: 

(2.2) y s=( c+ / g(r,y r ,z r )dr - \ z r dB r , 0<s<t, 

J s J s 

satisfies Y s G K , Vs G [Q,t], P-a.s.. 

(ii) y(t,z) G [0, T] x R nxd and\/y G R n such that d 2 K {-) is twice differentiable at y, 

(2.3) A{y-U K {y),g{t,U K {y),z)) < (D 2 d 2 K {y)z, z) + Cd 2 K (y), 



where C > is a constant which does not depend on (t,y,z), D 2 f{y) denotes the second order 
derivative of the function f. 
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For the following two BSDEs, 

(2.4) Vl = e+ [ 9 l (r, y l r , 4) dr- I z l r dB r , < s < t, i = 1, 2. 

J s J s 

Hu and Peng (2006) proved 

Proposition 2.2 Suppose that g 1 , g 2 satisfy (HI) ~(H3). Then the following are equivalent: 
(i) For any t G [0,T], V£\£ 2 G L 2 (T t ;R n ) such that ^ > £ 2 , the unique solutions (Y X ,Z X ) 
and (Y 2 ,Z 2 ) in S^(0,t;R n ) x 1^(0, t; R nxd ) to BSDEs (2.4) over time interval [0,t] satisfy: 

Y}>Y 2 , se[0,t], 

(n)Vt, V(y,y'), V(z,z>), 

n 

(2.5) -4(y-,g 1 (t,y + +y',z)-g 2 (t,y',z')) < 2 £ / yfe<0 |z fc - z' k \ 2 + C \y~ f , 

k=l 

where C > is a constant. 

We could not get an obvious intuition of a greater-than or less-than relation between g 1 and 
g 2 through (2.5). An explicit expression of (2.5) is given in the next section just like the well 
known one dimensional comparison theorem, see Theorem 2.2 in El Karoui, Peng, and Quenez 
(1997). 

3. EXPLICIT EXPRESSION FOR MULTIDIMENSIONAL COMPARISON 
THEOREM AND VIABILITY ON A RECTANGLE 

3.1. Comparison theorem 

This subsection presents a criterion equivalent to condition (2.5) given by Hu and Peng (2006). 
This criterion, of simpler form than Hu and Peng's condition has been exploited as a sufficient 
condition to prove the multidimensional comparison theorem by Zhou H. (1999) and Zhou S. 
(2004). By a short and forward calculus, we show that Hu and Peng's criterion coincides with 
condition in Zhou H. (1999). 

Theorem 3.1 Suppose that g , g 2 satisfy (HI) ~(H3). Then the following are equivalent: 

(i) For any t £ [0,T], V£\£ 2 G L 2 (T t ;R n ), such that > £ 2 , the unique solutions (Y x ,Z l ) 
and (Y 2 ,Z 2 ) in S^(0, t; R n ) x L%;(0, t; R nxd ) to BSDEs (2.4) over time interval [0,t] satisfy, 

Y}>Y 2 , se[0,t], 

(ii) For any k=l,2,...,n, Vi G [0,T], Vy' G R n , 

(3.1) gl(t,5 k y + y',z)>g 2 (t,y',z'). 

for any 5 k y G R n such that 5 k y > 0, (5 k y) k = 2 z, z' G R nxd and (z) k = {z') k . 

2 Here we follow notations in Hu and Peng (2006) for easy comparison. We denote by 5 y a constant belonging 
to R n associated with and (5 k y)k, the fcth component of 5 k y. 
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Proof. From Proposition 2.2, (i) is equivalent to 

n 

(3.2) -A(y-,g 1 (t,y + + y',z)-g 2 (t,y',z')) < 2 £ I Vk<0 \z k - z' k f + C \y~ | 2 . 

k=l 

Putting in (3.2) y = 5 k y + y k ■ e k , where y k £ R~ , S k y G R n , 5 k y > 0, {5 k y) k = 0, we get 

4y k (gl(t, d k y + y', z) - g 2 k (t, y' , z')) < 2\z k -z' k \ 2 + C \y k \ 2 . 

Taking z k = z' k , dividing Ay k and letting y k — > 0~, we deduce that g\{t, 5 k y + y', z) > g k (t, y', z'). 

Conversely, From assumption (H2) and condition (ii), we have Vi G [0, T], Vy' G R n , Vz,z' G 
R nx and z = (z[, z 2 , ...z' k _ 1 , z k , z k+1 , z' n ), 

fffc(*, $ k y + y', z ) - 9k(t, y', z') = g\(t, 5 k y + y', z) - gl(t, y', z") + g k (t, y' , z") - g k (t, y' , z) 

> - L \z k - 4| • 

Then we have for y = y k ■ e k , y k G R, y k < 0, 

4y fc (5fe(i,5 fc y + y / ,z)-^(t,y , ,z')) < ~^y k \z k - z' k \ < 2 \z k - z' k \ 2 + 2L 2 |y fc | 2 

(3-3) = 2/ yfc<0 |za; - 4| 2 + 2L > 2 \y~\ 2 1 

when getting the second inequality, we use 2ab < a 2 + b 2 . From (3.3) we deduce easily that 

n 

-4 (y-, g% y+ + y', z) - g 2 (t, y' , z')) < 2 £ I yk<0 \z k - z' k \ 2 + C \y~ | 2 . 

fe=l 

□ 

By Theorem 3.1 and its proof , we have the following corollary which has been given in Hu 
and Peng (2006). 

Corollary 3.1 Suppose that g 1 = g 2 = g satisfies (HI) ~(H3). Then the following are 
equivalent: 

(i) For any t G [0,T], V£\£ 2 G L 2 (F t ;R n ), such that ^ > £ 2 , the unique solutions (Y 1 ,^ 1 ) 
and (Y 2 ,Z 2 ) in Sj:(0,t;R n ) x L^O, t; R nxd ) to BSDEs (2.4) over time interval [0,t] satisfy: 

Y}>Yl 8€[0,t], 

(ii) For any k = 1,2, ...,n, for any t and y' , g k does not depend on {zj)-^ k and 
(3.4) g k (t,5 k y + y',z k ) > g k (t,y , z k ), 

for any 5 k y G R n such that 5 k y > 0, (5 k y) k = 0. 

Proof. In (3.1), let 5 k y = 0, z = (z±, z n ), z' = (0, z k , 0), and interchange positions 
of z' and z, we can deduce that for any t and y, g k does not depend on {zj).^ k . Then the 
equivalency of (i) and (ii) is a consequence of Theorem 3.1. □ 
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Remark 3.1 If condition (3.4) holds for g, we say that g(t, ■, z) is quasi-monotonously increas- 
ing. The "quasi-monotonously increasing" condition was first introduced for stochastic differen- 
tial equations by Mel'nikov (1983) and later cited by Geif3 and Manthey (1994), Ding and Wu 
(1998) to prove a comparison theorem for multidimensional stochastic differential equations. 
When condition (3.4) holds for g, by a direct calculation, the Lipschitz condition 

\9k (s,y,z k ) -g k (s,y',z' k ) \ < [i (\y - y'\ + \z k - z' k \) 

is equivalent to 

-^[^(yi-y'iy + \yk-y'k\ + \ z k- z' k \\ < g k (s,y,z k ) - g k (s,y',z k ) 

^ ^\^Z{yi-y'i) + + \yk-y'k\ + \ z k- z' k \ 

for each k = 1, . . . , n. 

3.2. Uniqueness theorem 

For two multidimensional BSDEs, if their generators coincide, of course they have the same 
solution; conversely if for any same terminal datum, they always have the same solution, do 
their generators coincide? We give an affirmative answer in Theorem 3.2. 

Proposition 3.1 Suppose that g satisfies (HI) ~(H3). Then the following are equivalent: 

(i) For any t G [0,T], V^ 1 ,^ 2 G L 2 {F t ]R n ) such that (f 1 - £ 2 ) G K, the unique solutions 
(y 1 ,^ 1 ) and (Y 2 ,Z 2 ) in S]r(0,t;R n ) x L^(0, t; R nxd ) to BSDEs (2.4) over time interval [0,t] 
satisfy: 

(Y}-Y 2 )^K, VsG [0,t], 

(ii) Mt, V{y,y'), V(z,z'), 

(3.5) 4 (y - Il K {y),g x {t, U K (y) + y' , z) - g 2 (t, y', z')) < (D 2 d 2 K (y)(z -z'),z- z') + Cd 2 K (y), 

where C > is a constant which does not depend on (t,y,z), D 2 f(y) denotes the second order 
derivative of the function f. 

Proof. Consider the following BSDE over time interval [0,t], 

(3.6) f s =l+ J^g (r, Y r , Z^ dr - J Z r dB r , 
where for y = (y 1 ,?/ 2 ), z = (z , z 2 ), 

~i \ I \, — 1 . — 2 — 1 ■ ~2\ 2/ ~2 ~2\ 2/ ~2 ~2\\ 

9{r,y,z) = (g {r,y +y ,z + z ) - g {r,y ,z ),g (r,y ,z )) . 
Set Y s = (3? - Y 2 ,Y 2 ), Z s = [Z] - Z 2 ,Z 2 ). Then (i) is equivalent to the following: 
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(iii) For any t G [0,T], V£ = such that G L 2 (T t ;K), the unique solutions (Y,.zQ 

in Sjr(0,t;R n ) x i£(0, t; R nxd ) to BSDE (3.6) over time interval [0,t] satisfy: 

Y, 1 G K, Vs G [0,t]. 

Applying Proposition 2.1 to BSDE (3.6) and the convex closed set K x R n , obviously (iii) is 
equivalent to (ii). □ 

As a consequence of Proposition 3.1, we have the following uniqueness result for generators of 
multidimensional BSDEs (2.4). 

Theorem 3.2 For any t G [0,T], V^ 1 = £ 2 G L 2 (T t ;R n ), 

Y s l = Y 2 , V*€[0,i], 
if and only if for any k = 1,2, n, V(t, y, 2) G [0, T] x R n x fl nxd , 

g\{t,y,z) = gl{t,y,z). 

Proof. Take K = {0}, by Proposition 3.1, It is seen that "Y 1 = Y 2 " is equivalent to 

n 

(3.7) 4 <y, ^(i, y>, z) - g 2 (t, y', z')) < 2 £ \z k - z' k f + C \y\ 2 . 

k=l 

Putting in (3.7) y = y k ■ e k , where y k G R, y k > 0, z = z' , we deduce that for any k = 1,2, n, 
y{t,y,z) G [0,T] x R n x R nxd , g\{t,y,z) < g 2 (t,y,z). Similarly the case ll y k < 0" leads to that 
V > g 2 "- Thus we obtain the desired result. □ 

Remark 3.2 Comparison Theorem 3.1 does not imply the above uniqueness result. If one ofY % 
is increasing in then it just implies that "£} = ^Y 1 = Y 2 44> g\{t,y,z k ) = g k (t,y, z k )", in 
which for any (t,y), g k does not depend on {zj)-^ k , due to the "quasi-monotonicity condition" 
(3.4). Note that in Theorem 3.2, g could depend on (zi, . . . , z n ). 

3.3. Viability on a rectangle 

Rectangles are special closed convex subset of R n . We first give explicit expressions for neces- 
sary and sufficient conditions for non-negative and non-positive solutions of BSDEs. 

Theorem 3.3 Suppose that g satisfies (HI) ~(H3). Then the following are equivalent: 

(i) For any t G [0,T], V£ G L 2 (F t ;R n ), such that £ > 0, (£ < 0), the unique solution (Y, Z) 
over time interval [0, t] satisfies, 

Y>0, (Y s <0),sG [0,t], 

(ii) For any k = 1,2, n, Vi G [0, T], 

(3.8) g k (t, 5 k y, 5 k z) > 0, (g k (t, -5 k y, 5 k z) < 0) 

for any 5 k y G R n such that 5 k y > 0, (5 k y) k = 0, and 5 k z G R nxd , (5 k z) k = 0. 



MULTIDIMENSIONAL RISK MEASURE 



9 



Proof. If (i) holds, then from Proposition 2.1, we have 

n 

(3.9) -4{y-,g(t,y+,z)) < 2 £ I Vk<0 \z k \ 2 + C 



y 



2 



k=l 



Take y = 5 k y + y fc • e fc in (3.9), y k G fl" 5 fc y G 12* <5 fc y > 0, {5 k y) k = 0, we get 

4y k g k (t,5 k y,z) < 2\z k \ 2 + C\y k \ 2 , 

dividing 4y k and letting z k = 0, y k — > 0~, we deduce that #&(£, 5 fc y, S k z) > 0. 
Conversely, if (3.8) holds, by assumption (H2), we have Vi E [0, T], 

5fc (t, 5 fc n, z) - g k (t, 5 k y, 5 k z) > -L \z k \ , 

thus g k (t,5 k y,z) > g k (t,5 k y,5 k z) - L\z k \ > -L\z k \. Then for y k G R, y k < 0, 

4y k {g k (t,S k y,z) < -4Ly k \z k \ < 2 \z k \ 2 + 2L 2 \y k \ 2 = 2I yk<0 \z k \ 2 + 2L 2 \y~ | 2 , 

Thus we can get easily that-4 (y~, g(t, y + , z)) < 2££ =1 I yk<0 \z k \ 2 + C\y~\ 2 . 

The result in the case "£ < 0" can be deduced from the preceding case "£ > 0" applied to the 
well chosen function g(t, y, z) := -g(t, -y, -z), V(t, y, z) G [0, T) x R n x R nxd . □ 

Minus C e R n on both sides of BSDE (2.2), we have 

(3.10) y s - C = £ - C + ! g(r,y r -C + C,z r )dr - ! z r dB r , < s < t, 

J s J s 

let Y s = Y s — C, £ = £ — C, y(r, y, z) = g(r, y + C, z), then apply Theorem 3.3 to the following 

(3.11) Vs=l+ g(r,y r ,Zr)dr - z r dB r , 0<s<t, 

J s J s 

we have 

Corollary 3.2 Suppose that g satisfies (HI) ~(H3). For any t G [0, T], 
(i) Assume C G R n . V£ G L 2 (J" t ; IT), aucft that £ > C, £ < <7j, tfien 

^>C, (y s <C7),VsG [0,t], 

if and only if for any k = 1,2,..., n, Vi G [0, T], 

(3.12) g k (t, 8 k y + C, 5 k z) > 0, (g k (t, -5 k y + C, 5 k z) < 0) 

/or any S k y G i? n suc/i i/mi 5 k y > 0, (5 k y) k = 0, and <5 fc z G i? nxd , (5 k z) k = 0. 
fwj Assume condition (3.4) and C\,Ci G i? n . V£ G [Ci,C2], 

y. g [Ci,c 2 ],v s g [o,t], 

if and only if for any k = 1,2,..., re, Vt G [0,T], 

(3.13) g k (t, 5 k y + C 1 ,z k = 0)>0> g k (t, -5 k y' + C 2 , 4 = 0), 

/or any 5 k y,6 k y' G 12 n swc/i that 5 k y,5 k y' G [0,C 2 - Ci], (£ fe y) fc = (5V)* = 0. □ 
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Remark 3.3 Assuming condition (3.4), let 5 k y = 5 k y' = 0, C\ = C2 in (3.13), we can deduce 
that V£ = C , C is a constant in R n , 

Y s = C,se [0,t], 
if and only if for any k = 1,2, n, Vi E [0, T], /or any t and C , 
(3.14) 9k (t,C,z k = 0) = 0. 

In fact the condition (3.4) is not necessary for the constancy (see Theorem J^-l(i) in Section 4)- 
4. MULTIDIMENSIONAL ^-EXPECTATIONS 

4.1. Basic properties 

A BSDE is in fact a dynamical mechanism of nonlinear expectation. It is natural to define 
multidimensional ^-expectation by multidimensional BSDEs. 

Definition 4.1 For each £ E L 2 (T t ;R n ), let (Y, Z) be the solution of BSDE (2.2), we define 
multidimensional g-expectation by £g[£] = Yq, multidimensional conditional g-expectation by 

^ < ./V Y t ,8i ill . 

The following properties hold for the multidimensional conditional ^-expectation. 

Theorem 4.1 Suppose that g satisfies (HI) ~(H3). For all t E [0,T], 

(i) Constancy. VC G R n , 

£ t g [C\T s ] = C,Vse[0,t], 
if and only if for any k = 1,2,..., n, Vi E [0, T], 

(4.1) 9k(t,C,z = 0) = 0, 

(ii) Monotonicity. For all £\£ 2 E L 2 (J r i ; R n ) such that > £ 2 , 

fjKV-] > etiflFs], vse [o,t\, 

if and only if for any k = 1,2, n, Vt E [0, T], Vy' E R n , g k does not depend on {zj)-^ k and 

g k (t,6 k y + y',z k ) > g k (t,y',z k ). 

for any 8 k y E R n such that 5 k y > 0, (5 k y) k = 0. 

(Hi) Positivity. V£ E L 2 (J r t',R n ) such that £ > 0, then 

el[t\T 8 ] >o,v s e [o,i], 

if and only if for any k = 1,2, ...,n, Vi E [0,T], 

5fc (i,5 fc y,(5 fe z) >0, 

/or any <5 fc y E i? n sucfr i/iai <5 fc y > 0, (5 k y) k = 0, and <5 fc .2 E i? nxd , ((5 fc z)fc = 0. 
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(iv) Time- consistency. For any < s < t < T , 

(v) Homogeneity. V£ G L 2 (T t ; R n ), Va G i2, 

S t g [ai\T s \=a£ t g [i\T s ]yse%ti 
if and only if for all k = 1, 2, n, Vt G [0, T] ; 

g k (t,ay,az) = ag k (t,y,z), 

for all (y,z) G R n x i? nxc! . 

Proof, (i) The sufficiency is obvious. We now prove the necessity. Let (Y s ) be the first part 
of solution (Y, Z) of 

(4.2) Y S = C + I g(r, Y r , Z r ) dr - I Z r dB r , < s < t, 

J S J S 

Minus C G R n on both sides of BSDE (4.2) and let F s = Y s - C, then (Y S ,Z S ) satisfies the 
following 

F s = + / g (r, F r + C, Z r ) dr- Z r dB r , 0<s<t. 

J S J S 

Take K = {0}, then by Proposition 2.1, V£ = 0, Y s = 0, Vs G [0, t] if and only if V(t, y, «), 

n 

(4.3) 4(y, 5 (t,C7,z)) < 2 £ |^| 2 + C \y\ 2 . 

k=l 

Let z = 0, (coordinate- wise), y = y k ■ e k in (4.3), G and — >• + , we get that for any 
k = l,2,...,n, Vt G [0,T], g k (t,C,0) < 0; if y fc < 0, and y k -> 0~, we have g k (t,C,0) > 0. Thus 
for any k = 1,2, ...,n, Vt G [0,T], g k (t,C,0) = 0. 

(ii) and (iii) are consequences of Corollary 3.1 and Theorem 3.3 respectively. The time consis- 
tency (iv) is due to the uniqueness of solution of BSDE (2.2). 

To prove (v), applying Theorem 3.2 to the following two BSDEs: 

Y s = £ + / —g (r,aY r ,aZ r ) dr — [ Z r dB r , a / 0, s G [0, t], 
Js a Js 

Y s = £ + ! g{r, Y r , Z r ) dr - I Z r dB r , s G [0, t], 

J S J S 

we obtain that for any f G L 2 (T t ;R n ), Y s = ^[o^|J" a ] = £%[€\T S ] = Y s , Vs G [0,t] if and only 
if g(t, ay, az) = ag(t, y, z). When a = 0, it is just the constancy. □ 

Remark 4.1 The constancy in Theorem 4-1 is for all constant in R n . In fact for any fixed 
C G R n , we still have £*[C|J-" S ] = C, Vs G [0, t] if and only if for any k = 1,2, ...,n, Vs G 
[0,T], g k (t,C,0) = 0. The proof is analogous. For instance, given C G R n , a g- expectation with 
generator g := r(y — C) + 0z has the property of constancy for only C. 



12 



YUHONG XU 



Remark 4.2 We do not impose in advance special conditions on multidimensional g- expectations 
in order to explore the general relations between g- expectations and their generating mechanism 
g, which is good for showing the essence of properties for multidimensional g- expectations. As 
we see above, multidimensional g- expectations are more complicated than the 1-dimensional one 
because they are related to different dimensions. 

4.2. Jensen's inequality 

Let cp = (cpx, . . . , ip n ) be a multiple convex function, i.e., p k ■ R — > R, k = 1,2, ...,n, is 
convex: Vx,y G R, VA G [0, 1], p k (Xx + (1 + Ay)) < \ipk(x) + (1 + A) <Pk(y)- This section shows 
that if Jensen's inequality holds for a multidimensional conditional ^-expectation, then the g- 
expectation consists in fact of n one dimensional ^-expectations with generator g satisfying 
positively homogeneous and negatively super-homogeneous conditions w.r.t. z. 

Definition 4.2 Jensen's inequality for multidimensional conditional g-expectation holds, if for 
any t G [0,T], for all £ G L 2 (T t ;R n ), convex function ip : R n -> R n such that <p(£) G L 2 (T t ;R n ), 
and for any k = 1,2, n, a.s., Vs G [0, t], 

(4.4) slM^m^M^M^)- 

Let be the payoff of a derivative at time T. Let £^.[6e] De its price at time 0. The above 
Jensen's inequality means that a convex transformation of the payoff at time T yields a 
higher price than the convex transformation of the initial price. So traders in financial markets 
should take care of this difference. For instance, consider a European call option with payoff 
max [(St — L) , 0], where St is the price of the underlying stock at time T, L is the strike price. 
If Jensen's inequality holds, then £ 9 [max [(St — L) , 0]] > max[£ 9 [St — L],0], which means that 
the price of the call option is greater than the nonnegative value obtained by a direct calculation 
of (S T - L). 

Jensen's inequality is in fact a comparison for solutions of two BSDEs. Let <pk(x) = x + a, x G 
R,a G R + . From the Jensen's inequality we get that + o| Ts] > £g k [Cfcl-^s] + a > [£fc|.F s ], 

V7c = 1,2,..., n, a.s., Vs G [0, t}. So it is natural to require the "quasi-monotonicity condition" 
(3.4) for the generator g. 

Remark 4.3 Jensen's inequality can be used to evaluate the nonlinear expectation operator Eg. 
A familiar example is £g k [\Ck\ [Fs] > l^fel-^sll; f or an V k = 1>2, ...,n, a.s., Vs G [0, T]. For 
a nonlinear expectation £ l g , we have to impose additional conditions on g to obtain Jensen's 
inequality. 

Theorem 4.2 Suppose that g satisfies (HI) ~(H3) and condition (3.4). Jensen's inequality for 
a multidimensional conditional g-expectation holds if and only if 

(i) for any k = 1, 2, ...,n, g^ does not depend on y and V(t, 2%) G [0, T] x R d , 
g k (t, Xz k ) = Xg k (t, z k ), VA > 0; g k (t, \z k ) > Xg k (t, z k ), VA < 0. 

Proof. If condition (i) holds, the multidimensional ^-expectation is in fact n one dimensional 
(/-expectation. Then by analogous proof of Jensen's inequality for one dimensional ^-expectation 
(see Jiang and Chen (2004)), we can obtain the multidimensional Jensen's inequality. 

We wish to emphasize some interesting aspects above the converse. Assume that Jensen's 
inequality for a multidimensional conditional ^-expectation holds. For any y G R n , take p y (x) = 



MULTIDIMENSIONAL RISK MEASURE 13 

x + y, x G R n . From the Jensen's inequality and the "monotonicity" of £*[-|J-" s ], we have 
Vs G [0, t], for any 8£ G L 2 (F t ] R%), 

efe + ss + ylFs] >^[^ s ] + y- 

Then respectively Y s := £g[£ + 5(, + y\Fs\ — V and Y s := <5*[^|J r s ] are solutions of 

(4.5) Y s = £ + 5£+ I g(r,Y r + y, Z r ) dr - f Z r dB r , < s < r. 

J S J S 

(4.6) Y s = £, + f g(r, Y r , Z r ) dr - I Z r dB r , < s < t, 

J s J s 

For BSDE (4.5) and BSDE (4.6), by Comparison Theorem 3.1, we have for any k = 1,2, ...,n, 
Vt G [0,T], for any y,y G iT, z fc G <5 fc y G iT and 5 k y > 0, (5 fc y) fc = 0, 

(4.7) 9k(t,y + y + 5 k y, z k ^j > g k (t, y, z k ) . 
Similarly taking (p- y (x) = x — y, then we can also get that 

(4.8) 9k{t,y~y + 5 k y, z k j > g k (t, y, z k ) . 

From (4.7) and (4.8), it is easy to deduce that g(t,y,z k ) = g(t,0,z k ). Therefore g does not 
depend on y. 

Take <p\(x) = Ax, A G i?\{0},x G R n . From Jensen's inequality and the monotonicity of 
£%[• \T S ), we have Vs G [0,t], VA G R\{0}, for any 5£ G L 2 {n, F t , P; R% ), 

£j[A£ + <5£|7y >A£j[e|7y. 

Then Y s := £*[A£ + 8£\T S ] and Y s := A£'[£|.F S ] are solutions of 



(4.9) 



Y s = A^ + 8£ + I g (r, Z r ) dr - I Z r dB r , < s < t. 

J S J S 



(4.10) Y s = A^ + J Xg (r, jZ r ^j dr - J Z r dB r , < s < r, 

For BSDE (4.9) and BSDE (4.10), by Comparison Theorem 3.1, we have for any k = 1,2, n, 
Vt G [0,T], VA G R\{0}, 

(4.11) g k (t,Xz k ) > Xg k (t,z k ). 
When A G (0, 1], we have 

(4.12) ^9k(t, z k ) > g k (t, jz k ). 
Thus VA G [l,oo), 

(4.13) Xg k (t,z k ) > g k (t,Xz k ). 
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Then combining (4.13) with (4.11), we get that VA G [1, oo), gk(t,Xzk) = Xgk(t,Zk). When 
A G (0, 1), the positive homogeneity holds too. 

Let Zk = in (4.11), we get gk(t,0) > \gk(t, 0), VA G i?\{0}. Therefore for any & = 1,2, ...,n, 
VtG[0,T], 

(4.14) 0k(t,O)>-Sfc(t,O), 
and 

(4.15) 5fc (t,0) > 2«fc(t,0). 

(4.14) and (4.15) lead to that for any k = 1,2, ...,n, Vt G [0,T], ^(t,0) = 0. □ 

Remark 4.4 The above theorem shows that if Jensen's inequality holds for a multidimensional 
g-risk measure (up to a minus sign, nonlinear expectation), then this risk measure is not recur- 
sive (g does not depend on y) in the generating mechanism g and different dimensions are not 
interacted with each other. The risk could be measured dimension by dimension respectively. 

We refer Jia and Peng (2010) for Jensen's inequality of g-convex function under one dimen- 
sional (/-expectation, g-convexity is a new notion of convex functions. They study this problem 
from a different point of view: for a given generator g, to characterize the function (p (not 
necessarily convex) such that the Jensen's inequality holds. For Jensen's inequality (4.4) and 
Jensen's inequality in Jiang and Chen (2004), we aim to characterize the generator g such that 
the Jensen's inequality holds for all convex functions. Note also that a convex function is not 
necessarily a ^-convex function. 

5. MULTIDIMENSIONAL DYNAMIC RISK MEASURES 

Since Artzner et al.'s seminal work on coherent risk measures, see Artzner, Delbaen, Eber, and 
Heath (1997, 1999), many authors are interested in the problem of how to measure the risk of a 
financial position as mentioned in the introduction. However they discover that the classical risk 
measure -VaR, could not satisfy some natural time-consistency. The notion of time-consistency 
was first proposed by Peng (2005) for nonlinear expectations (up to a sign, nonlinear risk mea- 
sures) in the framework of Brownian filtration. Recently a lot of papers concern time-consistent 
dynamic risk measures, see Bion-Nadal (2008a, 2009a), Cheridito, Delbaen, and Kupper (2006), 
Kloppel and Schweizer (2007) for characterizations of time consistency on a general filtered prob- 
ability space; see Peng (2005), Rosazza Gianin (2006), Jiang (2008) for dynamic g-risk measures 
introduced by BSDEs; see also Barrieu and El Karoui (2006), Jobert and Rogers (2008) and Ro- 
orda, Schumacher, and Engwerda (2005) for time consistent risk measures. This section defines a 
multidimensional risk measure by multidimensional ^-expectation and shows its properties and 
applications. 

Let (O, (J 7 t ) tg j T j , P) be a filtered probability space. Let Qt be the space of all financial 
positions with maturity T, i.e., the set of R n -v&\ued J^-measurable functions defined on f2. An 
element of Qt may be the net worth at maturity T of a financial contract, or the accounting 
value of a firm's equity, or the surplus of an insurance company. 

Definition 5.1 A multidimensional static risk measure is a functional p : Qt^R 71 such that 
(i) PiC 1 } — /°[£ 2 ]; ^/V^ 1 > £ 2 G Qt; (H) p[0] = 0. The quantity p[£] represents the riskiness of 
financial position £ G Qt- £, is acceptable if p[£] < 0, and unacceptable otherwise. 
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We call p Sj t '■ Qt^-L 2 {F s ;R n ), < s < t < T , a multidimensional dynamic risk measure if 
(Al) Monotonicity. Vs G [0,t], p s , t [^} < PsM% tfVf 1 > ? € ft; 
W atK] = -C, Ve € ft; 

(A5) Time- consistency. p s ,t[£,] = Ps,t[—pt,T[£]], V£ G ft 1 , VO < s < t < T; 
(A4) Regularity. l A p S)t [^l A } = lAPs,t[€\> V £ £ &t, VA G 
(A5J Normalization. p S) t[0] = 0. 

Condition (A5) means that "zero position, zero risk". Rosazza Gianin (2006) interpreted the 
time consistency as: computing the static risk measure p[£\ directly is equivalent to calculating 
p[£\ in two steps, i.e. evaluating first the riskiness pt,T[£] of £ at an intermediate date t and then 
quantifying at time the risk of — Pt,r[£] through p. The financial meaning of monotonicity is 
clear: The downside risk of a position is reduced if the payoff profile is increased. In a sense, the 
additional information available to the agent has to be fully used when assessing the riskiness of a 
payoff £, which means in particular that if we know that certain event A G T s is prevailing, then 
the riskiness of £ should depend on only what is really possible to happen, i.e. on the restriction 
of £ to A. This is characterized by (A4). The property (A4) is equivalent to: p s> t [£ 1a + £ 2 1a c ] = 
p,,t[?]lA + Ps,*[£ 2 ]1ac, VA G F s . Also, (A4) plus (A5) implies ft, t [?li] = Ups,t[^ v £ 6 ft, 
VA G J^. We refer to Detlefsen and Scandolo (2005) and Peng (2005) for the proofs. 

5.1. Multidimensional g-risk measure 
We call the generator g risk mechanism. In the sequence we consider Qt=L 2 {Ft, R n ). 

Definition 5.2 Suppose that the risk mechanism g satisfies (HI) ~(H3). For anyt G [0, T], the 
risk position £ G L 2 {Ft,R n ), we define p 9 Q t [£] = £g[— £], t [£] = £*[— ^.Fs], s G [0,t], as i/te muZ- 
tidimensional static g-risk measure and multidimensional dynamic g-risk measure respectively. 
Without ambiguity, we denote p 9 = Pq T , p 9 = p 9 tT . 

The financial meaning of pf[£] is: Let £ be the negative part of an uncertain income at time 
T, then — £ represents the potential loss at time T. Mathematically speaking, losses are viewed 
as negative gains. The ^-expectation is in fact a pricing operator (Peng (2006)). Therefore 
£]j[—£\Ft\ is the discount value at time t of potential loss — £. It is the risk for an investor. 
The fcth component p 9k [£, k ] of p\ [£] = (p 91 ^ 1 ], ■ ■ ■ ,P 9n [£, n ]) represents the risk of the kth sub- 
sidiary of a firm, pf 1 . . . , p 9n [£ n ] are related to each other through the generator (<?i, . . . g n ). 
The link between measures of risk and BSDEs is particularly interesting because it enhances 
interpretation and tractability of risk measures. Barrieu and El Karoui (2006) point out that 
the coefficient g of BSDEs can be interpreted as infinitesimal risk measure over a time interval 
[t,t + dt] as E[dYt\Ft] = —g(t,Yt, Zt)dt where Zt is the local volatility of the conditional risk 
measure, V(dYt\Ft) = \Zt\ 2 dt. Choosing carefully the coefficient g enables to generate dynamic 
g-risk measures that are locally compatible with the views and practice of different agents in 
the market. 

For multidimensional g-risk measures, additional condition has to be imposed to insure the 
monotonicity. Even for the convex and coherent risk measures in the next section, the quasi- 
monotonicity condition (3.4) is indispensable. It is easy to prove 

Theorem 5.1 Suppose that the risk mechanism g satisfies (HI) ~(H3) and (3.4) andg(-, 0, 0) = 
0. Then p 9 [-} and p 9 st [-] are multidimensional static risk measure and multidimensional dynamic 
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risk measure respectively. 

Proof. Theorem 4.1 yields that p 9 [£] is a multidimensional static risk measure and p% t [£] 
satisfies (Al) ~(A3) and (A5). We now prove (A4). Let (Y S ,Z S ), (Y S ,Z S ) be solutions of the 
following BSDEs 

(5.1) Y s = -^+f g(r,Y r ,Z r )dr - I Z r dB r ,0<s<t, 

J S J S 

(5.2) Y S = -£1 A + f g (r,Y r ,Z r ) dt - I Z r dB r , 0<s<t. 

J s J s 

Multiply 1a on both sides of the above two BSDEs, since V (y, z) G R n x R nxd , 

(5.3) lAg(r, y, z) = g(r, ylA, zl A ),\/A G F 8 , 
we get 

(5.4) Y s l A = -iU + [ g(r,Y r l A ,Z r l A )dr- [ Z r l A dB r , 0<s<t, 

J S J S 

(5.5) Y 8 1 A = -$1 A + ! g (r,Y r l A ,Z r l A ) dt - [ Z r l A dB r , 0<s<t. 

J s J s 

Comparing the above two BSDEs, by the uniqueness of solution of BSDE (— £1^, g, t), we obtain 
Y S 1 A = Y S 1 A , thus \aPs,AO-a] = lAPs,tVi], V^4 G T 8 . The proof is complete. □ 



5.2. Multidimensional convex and coherent risk measures 

Definition 5.3 A multidimensional dynamic g-risk measure pf t [-] is convex if for any t G 
[0,T], and for all risk positions ^ , £ 2 G L 2 (T t ;R n ), VA G [0,1], Vs G [0,t], p 9 s t [A^ 1 + (1 - A)£ 2 ] < 
A<je 1 ] + (1-AKJ^ 2 ]. 

Consider two strategies: one leads to the potential loss £ x , while another leads to the potential 
loss £ 2 . If one diversifies, spending only the fraction A of the resources on the first strategy 
and using the remaining part for the second alternative, the potential loss is A^ 1 + (1 — A)£ 2 . 
Thus, the axiom of convexity gives a precise meaning to the idea that diversification should not 
increase the risk. Monotonicity is a natural condition for one dimensional g-risk measures. Here 
we still need it for multidimensional g-risk measures. For a convex g-risk measure, we have 

Theorem 5.2 Suppose that the risk mechanism g satisfies (HI) ~(H3). For any t G [0, T], 
for all risk positions ^,^ 2 G L 2 {Ft,R n ), the multidimensional dynamic g-risk measure p 9 st [-] 
is nonincreasing and convex if and only ifit G [0, T], V(y%z l ) G R n x R nxd , i = 1,2, for all 
k = 1,2, ...,n, gt does not depend on {zj)-^ k and VA G [0,1], 

(5.6) g k (t,Xy 1 + (l-X)y 2 -d k y, z k ) < Xg k (t, y\ z\) + (1 - X)g k (t, y 2 , z 2 k ), 

for any 5 k y G R n such that 5 k y > 0, (5 k y) k = and z k = Xz\ + (1 - X)zf. 
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Proof. Let (Yf, Z\\ (Y 2 , Z 2 ), (Y S ,Z S ) be solutions of the following BSDEs 

Y} = -C 1 + f g (r, Y r \ Z\) dr - f Z l r dB r , < s < t, 

J s J s 

Y 2 = -£ 2 + f g ( r , Y T 2 , Z 2 ) dr- f Z 2 dB r , 0<s<t, 

J s J s 

Y s = -AC 1 - (1 - X)f + [ g {r,Y r ,Z r ) dr- I Z r dB r , 0<s<t. 

J s J s 

respectively. Set Y s = (AY/ + (1 - X)Y S 2 - Y s , Y s \ Y s 2 ), Z s = {\Z l s + (1 - X)Z 2 - Z s , Z\, Z 2 ). 
For y = (y\y 2 ,y 3 ) G (R n ) 3 , z = {z l ,z 2 ,z 3 ) G (R nxd ) 3 , A G [0,1], we define 

g(t, y, z) = (Xg(t, y 2 , z 2 ) + (1 - X)g(t, y 3 , z 3 ) - g(t, Xy 2 + (1 - A)y 3 - y\Xz 2 + (1 - X)z 3 - z 1 ), 
g{t,y 2 ,z 2 ), g(t,y 3 ,z 3 )), 

Then the convexity and monotonicity of p 9 t [-] is equivalent to the following 

(i) For any t G [0,T], V£ = (£\£ 2 ,£ 3 ) such that ^ > 0, the unique solutions (y s ,Z s ^ in 
S^(0, t; {R n ) 3 ) x L3r(0, t; (R nxd ) 3 ) to the following BSDE over time interval [0, t]: 

(5.7) Y a = l+ J g (r, Y r , Z^j du - j Z r dB r , 0<s<t. 

satisfies Y G -R™ x R n x R n . 

Applying Proposition 2.1 to the above BSDE and the set R+ x R n x R n , it is easy to see that 
(i) is equivalent to 

(ii) Vi, V(y,y'),V(z,z'), 

- (Ay-,Xg(t, y\ z 1 ) + (1 - X)g(t, y 2 , z 2 ) - g(t, Xy 1 + (1 - A)y 2 - y+, z)) 

n 

(5.8) < 2Y,Iy k <o\Xzl + (l-X)z 2 -z k \ 2 + C\y-\ 2 , 

k=l 

where C > is a constant which does not depend on (t, y, z). 
Similarly to Theorem 3.1, (ii) is equivalent to 

(iii) Vt G [0,T], for any k = 1,2, ...,n, VA G [0,1], 

g k (t, Xy 1 + (1 - A)y 2 - 5 k y, z) < Xg k (t, y\ z 1 ) + (1 - X)g k (t, y 2 , z 2 ), 

for any 5 k y G R n such that 5 k y > 0, (5 fc y) fc = and z, z 1 , z 2 G J R nX(l! , z fc = Xz\ + (1 - A)a£. 

Let A = 1, 5 k y = 0. we obtain that (?&(£, y 1 , z) < g k (t,y l , z 1 ), for any z,z x G i? nxc( such that 
Zk = -Zfc- Then it is easy to deduce that g k does not depend on {zj)-, k . Thus (5.6) holds true. □ 

Definition 5.4 A multidimensional dynamic g-risk measure p 9 t [-] is coherent if it satisfies 

(i) Positive homogeneity. V£ G L 2 (F t ; R n ), Va G R + , p^Jaf] =W>?,t[£],Vs G [0,t]; 

(ii) Subadditivity. For all risk positions £\£ 2 G L 2 (F t ;R n ), Vs G [O^f], p 9 st [^ + C } < P 9 S tl£ l ] + 

pLK 2 ]. 
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For a multidimensional coherent g-risk measure, we have 

Corollary 5.1 Suppose that the risk mechanism g satisfies (HI) ~(H3). The multidimen- 
sional dynamic risk measure p 9 st [-] is nonincreasing and coherent if and only if for all k = 

1.2, ...,n, Vi G [0,T], Vy G R n , g k does not depend on {zj)-, k , Va G R + , 

(5.9) g k (t,ay,az k ) = ag k (t,y,z k ), 
and V(y i ,4) G R n xR d ,i = 1,2, 

(5.10) g k {t,y l +y 2 -5 k y, z\ + z\) < g^y 1 , z\) + g k {t,y 2 , z\), 
for any 5 k y G R n such that 5 k y > 0, (5 k y) k = 0. 

Proof. Note that a coherent risk measure is equivalent to a risk measure which is convex and 
positively homogeneous and for a positively homogeneous function, it is subadditive if and only 
if it is convex. Thus this corollary holds. □ 

There are examples of dynamic risk measure on filtered probability spaces which are more gen- 
eral than those generated by a Brownian motion. See Bion-Nadal (2008a, 2009a) and Kloppel 
and Schweizer (2007). Even in the case of Brownian filtration, g-risk measures generated by 
Lipschitz drivers can not cover the whole. See the well known entropic risk measure introduced 
by Follmer and Schied (2002b) and Barrieu and El Karoui (2006) which is characterized by a 
quadratic BSDE. There are also coherent risk measure satisfying the axioms of positive homo- 
geneity and subadditivity but can not be expressed by a ^-expectation. Peng (2010) introduced a 
kind of sublinear expectation, called G-expectation (capital G), via a class of mutually singular 
probability measures and found some applications on coherent risk measure. See Bion-Nadal 
and Kervarec (2012), Nutz and Soner (2010) for further study on risk measuring under model 
uncertainty. 

5.3. Dual representation for multidimensional dynamic convex g-risk measure 

In this section x G R n means a row vector if without specification. Let * denote the trans- 
position of a vector or a matrix. For a matrix A mxm , we define exp{A} = YlnLo h.^™ • 
pf[-] be a convex g-risk measure defined in Definition 5.3. Let g be the generator such that 
for any k = l,2,...,n, g k does not depend on {zj)-^ k , quasi-increasingly convex in {y,z k ), i.e. 
satisfying (5.6) and (Y, Z) the solution of BSDE associated to the pair (£,g). We consider the 
Fenchel-Legendre transform of g k , k = 1,2, ...,n, 

G k (t,b,c)= sup [(y,b) + {z k ,c} - g k {t,y,z k )}, 

(y,z)eR n xR d 

Since g is convex and Lipschitz continuous, we have the duality relation 

(5.11) 9k(t,y,z k )= sup {(y,b) + (z k ,c) - G k {t,b,c)], 

where Dq ■= {(uj,t,b,c) e(lx[0,T]x R n x R d \G{uj,t,b,c) < +00} is the effective domain of 

G. The (u, t)-section of Dq is denoted by Dq . We denote by A the set of bounded pro- 
gressively measurable processes(/3, 7) with {[3 k ,^ k ),k = l,2,...,n valued in R n x R d such that 



MULTIDIMENSIONAL RISK MEASURE 



19 



E Jq \G k (t, (3t,jt)\ 2 dt < oo. The boundedness condition on A means that for any 7) G A, 
there exists a constant (dependent of (/3, 7)) such that 1/3^1 + 17^1 < C, t-a.e, P-a.s., k = 1,2, ...,n. 
Let us consider the family of linear generators 

g^(t,y,z k ) = (y,& fe ) + (z fc , 7t fe ) - G k (t, (3 k , 7* ), (/3, 7 ) G A 
Given (/3, 7) G .A, we denote by (Y^' 7 , Z^' 7 ) the solution to the following linear BSDE, 



(5.12) 



where we denote Z 
( (ZW) 
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P 



T 



Z s dB s , < t < T, 



( P n ■ 


■ P nl > 












\;Y = 






\ f3 ln ■ 











jdB 



( {l\dB) 
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Theorem 5.3 Suppose that the risk mechanism g satisfies (HI) ~(H3) and condition (3.4). The 
multidimensional dynamic convex risk measure p 9 t [-\ over [0, T] has the following representation: 
For a risk position £ G 1?{Tt'-, R n ), 

(5.13) P?[£]=ess sup {i? Q7 Kexp{/ T /3^}|J- t ]-Qf T (Q 7 )), 

08,7)6.4 I Jt ) 

where Q 1 is the probability measure with density process dLt = LtjtdBt, Lq = 1 and af T (Q' y ) = 



E& 



jf 1 G(s, p s , 7 S ) exp{ (3 u du}ds\ 7 t 



is called the penalty term. 



Proof. Observe from relation (5.11) and the "quasi-monotonicity condition" (3.4) that for 
all (0,7) G A, for any k = 1,2, ...,n, Vt G [0,T], V(y,z fc ) G i? n x 



(5.14) 



gk(t,5 k y + y,z k ) > g^ '(t,y, z k ). 



for any 5 k y G i? n such that 5 k y > 0, (S k y) k = 0. Then by Comparison Theorem 3.1, Yt > Y t ^, 
Vt G [0, T]. Furthermore, since G is convex with a linear growth condition on its effective 
domain, for each (t, u, y, z), the supremum in therelation (5.11) is attained at (b(t, y, z),c(t, y, z)) 
belonging to the subdifferential of —g, and so (|6|, |c|) is bounded by the Lipschitz constant of g. 
By a measurable selection theorem (see e.g. Appendix in Chapter III of Dellacherie and Meyer 
(1975)), since Y, Z are progressively measurable, we may find a pair of bounded progressively 
measurable processes (P,j) such that for any k = l,2,...,n, Vt G [0,T], V(y, z k ) G R n x R d , 

(5.15) 9k(t,y,z k )=g^(t,y,z k ) = (y,$) + (z k ,$) - G k (t,$,^). 

Then by the Uniqueness Theorem 3.2, these two solutions coincide, and so 

sup if' 7 > y/- 7 = Y t > sup y/ ,7 ,vt G [0,T]. 



08,7)6.4 



08,7)64 
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Then by solving the linear BSDE (5.12), we end the proof. □ 

The elements of A can be interpreted as possible probabilistic models, which are taken 
more or less seriously according to the size of the penalty cr(Q 7 ). Thus, the value pf[£\ is 
computed as the worst-case expectation taken over all models Q 1 and penalized by a^(Q 7 ). 
Note that each Q 7 is equivalent to the probability measure P with Radon-Nykodim derivative 
^jjp-p 1 * = exp{ — i Jq / y u 'y*du + y a dB u }. See Delbaen, Peng, and Rosazza Gianin (2010) for the 
representation of the penalty term of one dimensional dynamic concave utilities by the theory 
of one dimensional (7-expectation. The representation of coherent g-risk measure is a particular 
case of the above representation since it corresponds to the conjugacy function G = or +00. 

Corollary 5.2 Suppose that the risk mechanism g satisfies (HI) ~(H3) and condition (3.4). 
For the multidimensional dynamic coherent risk measure p 9 [-] over [0, T], for a risk position 
£ G L 2 ($7, Tt-, P'i R n ), the following representation holds: 



-£exp{^ j3 u du}\F t 

where Q 1 is the probability measure with density process dLt = LtjtdBt, Lq = 1. 



(5.16) /%[■] = ess sup E Q 



5.4. Cash additive risk measures 

Definition 5.5 A multidimensional dynamic g-risk measure p 9 st [-} is called cash additive risk 
measure if it satisfies the translation invariance: 

PlM + °\ = P"M -Cyte [0, t], VC G i?,V£ e L\T t - R n ). 

Theorem 5.4 Suppose that the risk mechanism g satisfies (HI) ~(H3). The multidimensional 
dynamic g-risk measure pf[-] is cash additive if and only if for all k = 1,2, ...,n, Vi G [0,r], 
\/z G R nxd , g does not depend on y. 

Proof. Set F t := p 9 ^ + C] + C, Y s := p 9 ^]. Then Y s and Y s are the solutions of 
(5.17) Y s = -i+ f g (r,Y r - C,Z r ) dr - f Z r dB r , 0<s<t. 

J S J S 



(5.18) Y s = -£+ g(r,Y r ,Z r )dr - Z r dB r ,0<s<t, 

J s J s 

For BSDE (5.17) and BSDE (5.18), by the Uniqueness Theorem 3.2, we obtain the necessary 
and sufficient condition: for any k = 1,2, ...,n, Vt G [0, T], V(y,z) G R n x R nxd , 

(5.19) g{t,y-C,z)=g{t,y,z). 



Due to the arbitrariness of (y, z, C), we deduce that g does not depend on y. □ 

The above cash additivity also holds conditionally: G i 2 (J r t; R n ), V77 G L 2 (J r s ; R n ), p 9 s t [£ + 
v] = Psti£] — ^, Vs G [0, t] iff g does not depend on y. The cash additivity property is also called 
translation invariance. While the convexity and the monotonicity axioms have been largely ac- 
cepted by academics and practitioners, the cash additive axiom has been criticized from an 
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economic viewpoint. A basic reason is that while regulators and financial institutions determine 
and collect today the reserve amounts to cover future risky positions, the cash additivity re- 
quires that risky positions and reserve amounts are expressed in the same numeraire. This is 
a stringent requirement that limits the applicability of cash additive risk measures. Moreover, 
the cash additive axiom does not allow to account for stochastic discount factor. El Karoui and 
Ravanelli (2009) relaxed the cash additivity axiom to cash subadditivity. By Theorem 5.4, a 
multidimensional (/-risk measure satisfying the cash additive axiom and the monotonicity axiom 
is in fact n one dimensional g-risk measures put together. For a g-risk measure with genera- 
tor g := —rty+ g(t,z), where g(t, 0) = 0, the corresponding g-risk measure pfy[-] satisfies the 
following translation invariance axiom: Vi € [0, T], VC E R n ,(, £ 1?{Tt\ R n ), 



(5.20) 




5.5. Model uncertainty and robust measurement of risks 

When making a decision, different traders may use different probability measures. When mea- 
suring the risk of a financial position, we may face several models to choose. In these cases, we 
are interested in knowing the minimum and maximum of a family of risk measures. 

Assume that there is a family of multidimensional g-risk measures p 9 t with parameters (g a , £ a ) 
satisfying (HI) ~(H3) and condition (3.4). We will show that, under some mild conditions, the 
supremum ess suppf t is also a g-risk measure with parameters (ess supg a , ess sup^ Q ). 

a a a 

Theorem 5.5 Suppose that there is an a such thaty(y, z), P-a.s., t-a.e. 
(5.21) g(t, y, z) = ess supg a (t, y, z) = g a (t, y, z), 

a 



(5.22) f = ess sup£ Q = C- 

a 

Let p 9 st [£] be the solution of BSDE (g,£)- Then 

p% t [£] = ess suppg m = P 9 Z [r] ,W € [0, T], P-a.s.. 

a 

Proof. Since (g a ,£ a ) and (<?,£) satisfy the requirement of Comparison Theorem, for any a, 
we have p 9 st [£] > p 9 st and hence p 9 st [£] > ess suppf t Vi 6 [0, T], P-a.s.. 

a ' 

On the other hand, since (<?,£) is attainable at some a, p 9 st [£] and p 9 st are both solutions 
of BSDE (g a ,£ a ). By the Uniqueness Theorem we get that 

ess suppg m < Pit [£] = Ptt [e] < ess suppg • 

a a 

Hence we obtain the desired result. □ 

This theorem still holds when "sup" is replaced by "inf " . It is a generalization of dual repre- 
sentation for convex and coherent g-i\sk measures in which {g a } are linear generators. Thus we 
conclude that 
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• Convex and coherent g-risk measures are robust risk measures; 

• Model uncertainty could lead to the convexity of risk measures. 

Of course, model uncertainty we consider here is under one same probabilistic framework (or 
different probability measures which are equivalent to the reference probability measure). Peng 
(2010) considered coherent risk measures with different probabilistic models which are mutually 
singular in the language of G-expectation. See also Bion-Nadal and Kervarec (2012) for convex 
and non-dominated risk measures. 

5.6. Motivations and Applications 

Multidimensional BSDEs have been found several practical applications, among others, see 
mean- variance hedging in Mania and Tevzadze (2003), nonzero game in El Karoui and Hamadene 
(2003), switching problem in Hu and Tang (2008), Hamadene and Zhang (2010) and references 
therein. In this section we shall show how our multidimensional approach can be applied to 
measure the insolvency risk of interacted business units. 

5.6.1. Measuring the insolvency risk of a firm with interacted subsidiaries 

One way to measure the insolvency risk of a firm is to use the probability of insolvency itself. 
It is to determine the minimal amount of capital to add to the firm such that the probability 
of insolvency is less than or equal to some given probability (3. That is to find the minimal a 
such that P {X + a(l + r) < 0} < (3, where X is the risky asset hold in time T for the firm, 
r is the interest rate between and T. The solution to the above problem is called Value at 
Risk (VaR). It is shown that VaR is not a time consistent risk measure by Artzner, Delbaen, 
Eber, and Heath (1999). They studied various common risk measures including VaR. However 
they omitted an intuitive risk measure that one could construct: the cost of an insurance policy 
protecting the firm against insolvency. Jarrow (2002) defined the premium of a put option on 
the firm as a measure of insolvency risk. We recall Jarrow's idea detailedly. A firm is in fact 
a portfolio of assets and liabilities. The firm's net value is the value of the firm's assets minus 
liabilities. A measure of insolvency risk is the cost of buying a put option written on the firm' 
net value with zero strike price. Such a put option could guarantee that the firm's net value will 
always be greater than or equal to zero. 

Let us be more precise. Let (O, P) be a probability space. Let Qt denote the space of all 
financial positions with maturity T, i.e., the set of i? n -valued functions defined on f2. Let r 
denote the interest rate between time and T. A firm is defined by its balance sheet at time T. 
The balance sheet consists of assets A £ Qt-, liabilities C 6 Qt, and equity £ £ Qt such that the 
following accounting identities hold: 

(5.23) £ = A-C, if .4 - £ > and £ = if .4 - £ < 0. 

We call X = A — C the firm's net value. Note that it maybe negative while £ cannot due to its 
limited liability status. We say the firm is insolvent if X is strictly negative, or else the firm is 
solvent. Capital relates to the amount of the firm's net value invested in the riskless asset (e.g. 
cash). It is a different notion from equity. Let X be the risky assets minus liabilities and a be 
the firm's capital invested at time 0. Then the firm's net value X = X + a(l + r). We want to 
quantify the insolvency risk of the firm's net value X = X + a(l + r). As illustrated above, VaR 
is a such measure. Now we consider a put option written on the firm's net value X with strike 
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price zero. The terminal condition is 

(5.24) max[-(X + a(l + r)),0]. 

Under the condition of arbitrage-free, Jarrow (2002) defined the "Put Premium Risk Measure" 
as 

E I max 

(5.25) p[X] = — 

It measures the economic cost of avoiding insolvency. This is a single period (static) model. It 
can be generalized to the dynamic situation. 
Recall that the classical Black-Scholes formula corresponds to the following BSDE: 

(5.26) dY t = r t Y t + t Z t dt + Z t dB t , Y T = £, t€ [0, T], 

where (rt) is the instantaneous interest rate, (0t) is called a risk premium and (Z%) is viewed 
as the portfolio of risk investment. The meaning of (5.26) is that the change of the value of an 
option equals to the interest income plus the return of risk investment and plus a random noise. 
See El Karoui, Peng, and Quenez (1997) for a derivation. We then can define a dynamic g-risk 
measure pf T [£] with £ = max [— (X + a(l + r)) , 0] and g (t, y, z) = — (r t y + Otz). Generally in 
a financial market (may or may not be complete), we can define a more general g-risk measure 
as the solution of the following BSDE: 

(5.27) dY t = -g(t,Y t ,Z t )dt + Z t dB t , Yt = max [— (X + a(l + r)) , 0] , t e [0,T], 

The above put option provides insurance against the firm's insolvency. When the firm is insolvent, 
it pays off and covers the loss X + a(l + r) < 0. See El Karoui, Peng, and Quenez (1997) for an 
abundance of examples corresponding to various g. 

The above put premium risk measure can be summarized to more involved cases. Economic 
agents such as consumers, producers, banks, and stockholders, constantly interact with each 
other in different ways and for different purposes. They may compete, imitate and communicate 
with each other or suspect strategies of their counterparties and somehow out of these indi- 
vidual interactions a certain coherence at the aggregate level develops (Kirman (1997)). The 
first contribution in which the problem of stochastic interaction was explicitly treated is that of 
Follmer (1974). He showed that if the characteristics of agents, for example their preferences are 
random but dependent on those of others, the effect of large numbers of agents is not enough to 
eliminate uncertainty at the aggregate level. 

Indeed economic agents are continuously learning and changing their strategies. One is inter- 
ested in knowing how to measure risks they encounter, in particular the insolvency risk for a 
firm with a number of interacted subsidiaries. See the following examples. 

Example 5.1 Consider a firm with two subsidiaries competing against each other for finite 
resource (Lt) (e.g. cash income of the firm) at each time t to hedge their possible net values £ 
and £ 2 at time T respectively. Let (rj,r^) and (Q\,d\) be their risk exposure to the spot interest 
and risky investment. Let (Y^,Y^) and (Z\, Z%) be their hedging processes and portfolios of risk 
investment. The first agent suspects his counterparty may get a^Y^ from the finite resource (Lt) 
with proportional coefficient < < 1. Thus he gets (Lt — alYt) + . And the second agent 
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suspects his counterparty, the first agent may get a\Y^ from the finite resource (Lt) with propor- 
tional coefficient < a\ < 1. Thus he gets (Lt — a\Y^-) + . To reduce their insolvency risks, the 
two subsidiaries can decide whether or not to buy a put option with payoff max [-f ,0], i=l,2 
at time T. Then a natural problem comes up: how to estimate this put option? Using J arrow's 
formula (5.25) or the classical Black-Scholes equation (5.26) ? Obviously different pricing for- 
mulas or pricing under various markets ( complete or incomplete ) yield different current values. 
The current value is the capital a subsidiary needs to protect herself against insolvency. An ap- 
propriate and efficient way is that the two subsidiaries estimate the put options by themselves 
according their own flows of cash payments / incomes. That is, they spend the required capitals 
(equal to prices of the put options) on investment such that the risk capitals hedge the insolvency 
risks at time T. Therefore in a complete financial market the pricing processes of the two put 
options should satisfy 

(5.28) dY? = (r^ 1 + (L t - a 2 t Y t 2 ) + + Q\z\) dt + Z\dB u Y^ = max [-?, 0] , t € [0, T], 

(5.29) dY t 2 = (r 2 t Y 2 + (L t - a\Y^) + + 6 2 Z 2 ) dt + Z 2 dB t , Y$ = max [-£ 2 , 0] , t € [0, T]. 

Then (pi p 2 [£ 2 ]):=(Y t , Y 2 ) is a 2-dimensional risk measure but not sublinear. Y^ and Y 2 
are interconnected in their generators because of their reciprocity of benefits. The firm should 
allocate corresponding capitals (equal to p\, p 2 resp.) at time t to her subsidiaries to protect them 
against insolvency. 

See also the following simple example. 

Example 5.2 Consider a multi-business firm with two business units holding stocks of each 
other for sharing their Profit&zLoss. The first unit holds 20% stocks of the second one and the 
second unit holds 10% stocks of the first one. Thus their insolvency risks can be measured by 

dY t l = [r t (0.9^ + Q.2Y 2 ) + Q t Z\\ dt + Z\dB u Y^ = max 0] , 
dY 2 = [r t (O.l^ 1 + 0.8lf) + 6 t Z 2 } dt + Z 2 dB t , Y$ = max [-£ 2 , 0] . 

It is interesting that Stadje (2010) extended various risk measures in discrete time to the 
solutions of BSDEs, including VaR, semi-deviation risk measure, Average VaR and Gini risk 
measure. Thus, BSDEs provide an abundance of dynamic valuations (resp. risk measure) in 
continuous time. For instance, the entropic risk measure is the most famous example of a convex 
risk measure, introduced by Follmer and Schied (2002b). Barrieu and El Karoui (2006) pointed 
out that the classical dynamic entropic risk measure under a Brownian filtration corresponds to 

1 1 2 

a quadratic BSDE with generator g = \z\ , where 7 > is the risk tolerance coefficient. 

Along with the rapid development of the financial market, risk relationships will be more 
and more complicated. For instance, a set of related options, a group of players in a game, 
companies in the same industry, due to countless ties among them, it is more applicable to 
consider multidimensional risk measures to control their risks. 
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5.6.2. Acceptance sets and insurance g-risk measures 

From an insolvency point of view, an acceptance set identifies those risk positions that are 
acceptable such that the firm is solvent. Any risk measure p : Qj t-^R n induces an acceptance 
set 

A P :={teg T \ p[£] <0}. 
Conversely, for a given acceptance set A, the associated risk measure can be defined as 

p[t] := inf {a G R n \ £ + a G A] . 

Within the static framework of risk measures, Artzner, Delbaen, Eber, and Heath (1999) pre- 
sented four axioms for an acceptance set should satisfy. And they summarized the relations 
between a coherent risk measure and its acceptance set. In this section we put emphasis on the 
following definition: 

(5.30) pf iT [£] := ess inf {r? G L 2 {F t] R n )\ p 9 ^ + rj\ < o} , £ G L 2 (F T ;R n ), 

where p 9 tT is a g-risk measure. The measure p^ T is the minimal capital required to protect the 
firm against insolvency at time t. We call it "Insurance Risk Measure" as in Jarrow (2002). Of 
course if we define the following acceptance set 

A% := {^L 2 (T T ;R n )\p 9 T [(]<o}, 

then p 9 T can be also written as 

% >T [Z] ^essinf^GLV^i?")! £ + r, e A f pg } , £ E L 2 (F T ; R n ). 

If p 9 tT is convex in the sense of Definition 5.3, it is easy to prove that A* 9 is convex, i.e. 
V^ 1 ,^ 2 G A* 9 , VA G [0, 1], we have (A^ 1 + (1 - A)£ 2 ) G A* 9 . We are interested in the relation 
between p g tT and J? tT - Obviously we have 

Proposition 5.1 Let p? T be a multidimensional dynamic g-risk measure. Then 

(i) p^ T is always translation invariant. 

(ii) jof j, = p 9 tT if and only if g does not depend on y. 

(Hi) f^ T > p 9 T if and only if p 9 tT satisfies the axiom of cash subadditivity: 

(5.31) plAt + v] > P g t, T lZ]-V,Vv G L 2 {F t ;R n )^ G L 2 {F T -R n ). 

The static risk measure p^ T induced by p^ T performs even "better" than the original one be- 
cause it is always translation invariant and constant-keeping. It extends the notion of "Insurance 
Risk Measure" in Jarrow (2002) to generally nonlinear case. 

It seems to be more reasonable to establish the axiom of cash subadditivity as 

(5.32) P B t,T^ + r l }>p 9 t)T ^]- V ,y V EL 2 (T t ;Rl),^ G L 2 {F T - R n ). 

An example in the following (Example 5.3) shows that the typical risk-neutral BSDE verifies 

(5.32) but not (5.31). See Barrieu and El Karoui (2006) for other motivating examples. Now we 
consider the following (nonnegative) insurance g-risk measure: 

(5.33) -essinfj^GL 2 ^;^)! p 9 ^ + rj] < o} , £ € L 2 (F T ; R n ), 
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Theorem 5.6 Let g satisfy (HI) ~(H3) and (3.4) and <?(-,0, 0) = 0. Assume also that p\ T 
satisfies the cash subadditivity (5.32). Then (i) p t = ; V£ > 0; (ii) p^^[-] is nonincreasing; 
(Hi) > % T [£], V£ G L 2 (.F T ;i?«); (iv) pt^ + v] > P$[£]-V,VV G L\F t ;R%)M G 

The proofs are trivial. It is interesting that may be not time consistent. 

Example 5.3 Consider the 1- dimensional risk-neutral BSDE (£,g,T) with g{y) = —ry, r G 
R + , y G R. The induced g-risk measure p 9 tT [£\ = e~ r<yT ~^ ■ E[—^\Ft] satisfies the cash subaddi- 
tivity (5.32). It is easy to deduce that p^j! [£] = (E[— £| J~t]) + and it satisfies Theorem 5.6. For 
£ = — Bt and t G (0,T], we /iaue 

pjf [-^[-Br]] = p# [-Bt] = E[B+] > = (E[B T )) + = p$ [-Br] . 

TTius /5^j? violates the axiom of time consistency for some random variables. 

Risk management is different from the firm's and the regulator's perspective. Prom the point 
of view of supervision, regulators prefer p. £ to p 9 tT while the firms may just choose p? T to 
prepare their risk capital. 

5.6.3. Optimal risk sharing for ^-tolerant g-risk measures 

The topic of optimal risk sharing has been studied early by Borch (1962) in the insurance 
and reinsurance context and recently by Barrieu and El Karoui (2005, 2006) for cash additive 
7-tolerant convex risk measures. For non-homogeneous convex risk measures, the impact of the 
size of the position is not linear. It seems natural to consider the relation between risk tolerance 
and the perception of the size of positions. The risk tolerance coefficient 7 describes how agent 
penalize compared with the root risk measure p. Precisely a 7-tolerant risk measure is defined 
as 

7 

A typical example is the entropic risk measure eJ T [£] = 7ln£'[exp(— -)|J r t]- 

The optimal risk sharing between two economic parties, denoted by A and B means that: A 
want to issue a financial product r\ G L 2 (Ft) and sell it to B for price nt at time t. A tries to 
minimize her risk 

ess inf {ptr[£ — T]] — fit} subject to pfrhl] + Tit < 0. 

The price of rj is determined by B as irt = —pfj<[v] - Hence B should have an interest in doing the 
transaction. At least the product rj does not worsen her risk. Then the minimization problem 
becomes 

ess inf {pf iT [£ - rj\ + p^ T [rf\} , 

which is called inf-convolution of two functionals and denoted by p Up B [£]. 

Now we consider two groups A and B, measuring their risks by a multidimensional g-risk 
measure p 9 with different risk tolerant coefficients 7Ai7b G R + , i.e., their risks are measured by 
Ptjr '■= lAPt iT [^] and PtJ '■= lBPi, T i^] respectively. 
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Theorem 5.7 Let p 9 T be a multidimensional convex g-risk measure with g satisfying (HI) 
~(H3) and (3.4) and g(-,0,0) = 0. T/ien V7a,7b G R+ , 

(i) p g fA\3p9iB is a (7 a + 7b) -tolerant g-risk measure, i.e., it is the unique solution of the 
following BSDE 

(5.34) Y t = -e + ^ (7A + 7b) 9 (s, Y ' , Zs ) - C Z s dB s , < t < T. 

(ii) (p 9 "')^ =00 is a goo-risk measure with generator g^ = lim jg [t, ) . And for different 
g,g satisfying the above assumptions, we have 

p9°°n^ = { P 9 a P 9 ) oo . 

Proof, (i) By the convexity of p 9 tT and p 9 T [0] = 0, it is easy to check that p 9l T [£,} is non- 
increasing in 7 and pf T (7,£) := pp T [C] is sublinear (positive-homogeneous and subadditive) in 
(7,0- Thus for any n G L 2 (T T ; R n ), 

(5.35) ( 7A + 7b) pI.t L f ^ J < 7Apf, r [^— — ] + 7Bpf, r [^-], 

' (7A + 7Bj 1A 'IB 

and consequently (7A + 7b) P?,t[ ( 7a | 7b ) ] ^ ess ™ f {wj r [^?] + Tb^tI^]}- When taking 
f? = ( 7 J+ 7B ) , inequality (5.35) becomes an equality. Hence (7A + 7b) P g [ ( 7A | 7a) ] = P 97 - 4 Dp 97 ^ . 
Obviously (7a + 7b) pf t[ ( 7a + 7b ) ] ^ s un il ue solution of BSDE (5.34). Furthermore, rf = 



is the optimal transfer. 



(7A+7s) 

(ii) Since g is Lipschitz (thus dominated by positive- homogeneous functions, see Remark 3.1) 

and both p 9 ~* and g 7 := 75 (t, 9 , ^ are nonincreasing in 7, (700 and p 9 °° are well defined. Then by 

Theorem 5.5, p 9ca is the solution of BSDE with generator g^. Since p 9 " 1 and p 9 ^ are nonincreasing 
in 7, by a direct calculation, we have 

— £ — Tj — Tj 

p 9o °np 9o ° = ess inf ess inf{7p? r [ ] + 7P? r [— 1} 

rj 7 7 7 

= ess inf ess inf{7(p? T [— — 77] + pfyt 7 ?])} 

7 TJ ' 7 

= ess infp^Dp^ = (p^Up 9 )^ . 

□ 

The financial meaning of (i) is obvious: the behavior of risk sharing increases A's risk tolerance 
capacity which equals to the sum of the risk tolerance of both groups. 



5.7. Other ways to obtain g-risk measures 
There is another way to define risk measure via ^-expectation. 

Definition 5.6 Suppose that the risk mechanism g satisfies (HI) ~(H3) and (3.4) andg(-, 0, 0) = 
0. For any t G [0,T], the risk position £ G L 2 (T t ;R ), we define p 9 [£] = -£*[£], p 9 st [£,] = 
-^[^]^G[0,i] 
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Obviously p 9 [£,] and are multidimensional static risk measure and multidimensional 

dynamic risk measure respectively. However they have different meanings compared to y-risk 
measures in Definition 5.2. For a future risk position £ with maturity time t, if we give the price 
of £ via ^-expectation and view £g[£] as the estimation value, then the risk measure p~ 9 [£] means 
to measure the risk of £ via the opposite value of its estimation £1 [£] . A natural question comes 
up: for a specific risk position £, we use which mechanism to measure its risk, or equivalently 
which mechanism is more robust. 

5.7.1. Multidimensional g-pricing mechanisms 

In an option market, the ask-bid pricing mechanism is operated through the system of market 
makers. The ask price and the bid price often do not coincide with each other (see Jouini (2000), 
Bion-Nadal (2009b)). If £j[£\Tt] is the ask price at time t of a derivative £ with maturity 
T, then the bid price must be —£j[—^\Tt] and we have generally fj^l^t] > — £j[~ £| ^t]- A 
^-expectation is called a g-pricing mechanism here. We have 

Theorem 5.8 Let (HI) ~(H3) hold for g. For any t G [0,7], V£ G L 2 (T t ;R n ), £*[£|.F S ] > 
Sll-^Ts], Vs G [0,t], if and only if for any k = 1,2, ...,n, Vi G [0,7], Vy' G R n , 

(5.36) g k (t, 5 k y + y', z) > -g k (t, -y , -z), 

for any 5 k y G R n such that 6 k y > 0, (5 k y) k = 0, and z, z' G R nxd with (z) k = (z') k . 

Proof. Let (Sg[— £1^]) denote the solution of the following BSDE: 



(5.37) Y s = -Z + I g (r, Y r , Z r ) dr - f Z r dB r , 0<s<t. 

J S J S 

Then (-£*[-£ I-FJ) solves the following BSDE: 

(5.38) Y s = £ + j ~g (r, -Y r , -Z r ^j dr - j Z r dB r , 0<s<t. 

Applying Theorem 3.1 to the above BSDE (5.38) and BSDE (2.2), it is deduced that for all 
derivative £, > — £"*[— C| Fs], Vs G [0,t] if and only if for any k = 1,2, ...,n, Vt G [0,7], 

Vy' G R n , g k {t, 5 k y + y', z) > -g k (t, -y', -z'), for any 5 k y G R n such that 5 k y > 0, (5 k y) k = 0, 
and z,z' G R nxd with (z) fc = (zOfe- □ 

An axiomatic approach in continuous time asset pricing is developed in Jouini (2000). The 
axiomatization of dynamic pricing procedure in the setting of Brownian filtration is presented in 
Peng (2005, 2006) in the language of y-expectation. Pricing via dynamic convex risk measures 
has also been studied by Jobert and Rogers (2008). Bion-Nadal (2009b) introduced a notion of 
time consistent dynamic pricing procedure on a filtered probability space which assigns to every 
essentially bounded financial position a dynamic ask price and a bid price, taking transaction 
costs and liquidity risk into account. See also Bion-Nadal (2008b). Characterization of one 
dimensional g-pricing mechanisms by its generating function was studied sufficiently in Peng 
(2006). For a risk mechanism g, if inequality (5.36) holds, then we have > pf[£], V£ G 

L 2 (J : T',R n )- Thus usually we utilize /?f[£] rather than pf[£] to measure the risk of a financial 
position £ in an option market. 
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5.7.2. Utility-based g-risk measure 

Let us consider a small agent who can consume between time and time T. Let q be the 
(positive) consumption rate at time t. We assume that there exists a terminal reward £ at time 
T. The utility at time t is a function of the instantaneous consumption rate c± and of the future 
utility (corresponding to the future consumption). In fact, the recursive utility is assumed to 
satisfy the following BSDE in form of conditional expectation (Duffle and Epstein ( 1992a, b)), 

ut = E[l- + J g(s,c s ,u s )ds\J r t }. 

We consider a more general class of recursive utilities defined by Peng's BSDE (El Karoui, Peng, 
and Quenez (1997)), 

(5.39) u t = £+ / g (s,c s ,ul . . .u™,z s ) ds - z s dB s . 

Jt Jt 

In general, utilities must satisfy the following classical properties: (a) Monotonicity with re- 
spect to the terminal value and to the consumption, (b) Concavity with respect to the con- 
sumption, (c) Time consistency: this means that, for any two consumption processes c x and c 2 
and any time t, if c 1 and c 2 are identical up to time t and if the continuation of c 1 is preferred 
to the continuation of c 2 at time t, then c 1 is preferred to c 2 at time 0. Certainly by Theorem 
3.1 and Theorem 5.2, if g is concave w.r.t (c,y,z), quasi-monotonously increasing w.r.t (c,y), 
and for any (c, y), the generator g^ does not depend on (zj)j^ k , then the above properties are 
satisfied. From a financial point of view, the generator g(t,Ut,Zt) represents the instantaneous 
utility at time t of consumption rate cj. The dependence of on u n can be interpreted 

as a nonzero-sum game problem, where the players' utilities affect each other and consequently 
the generators are interconnected. 

Let Ut be a function of utility, a utility-based risk measure is defined by pt = —u% (see Kulikov 
(2007) for example). Hence we define the utility-based g-risk measure pf[-} by the solution u\ 
of BSDE (5.39). Note that unlike p\ [•], the convexity (subadditivity) of risk measure pf [•] is 
not consistent with its generator g, i.e.: if g is concave (superadditive), then pf [•] is convex 
(resp. subadditive). Thus for investigating convex utility-based risk measure, one usually studies 
concave utility (Kulikov (2007)). 

Observe that by Theorem 5.4, when g does not depend on y, the definition ~p g tT 
coincides with the following 

(5.40) pf r [£] := essinf { V € L 2 (T t ; R n )\£j + V \T t ] > 0} . 

So generally a risk measure based on a recursive utility is not consistent with definition (5.40). 
Based on some experimental evidence, Brachinger and Weber (1997) argued that in general a 
risk measure should not be considered, up to a sign, an utility functional. Let u be a concave, 
strictly increasing utility function with u(0) = 0. The risk measure p u is defined as p u [£] := 
inf {a G R\ Ep [u (£ + a)} > 0}. Frittelli and Scandolo (2006) also gave a partial support which 
shows that the preference structure defined by Epu does not coincide with the reverse ordering 
provided by p u , unless u is linear or exponential. Time consistency is a crucial property for 
utility-based risk measures. Cheridito and Kupper (2009) gave some equivalent conditions for 
time consistency of a dynamic utility function. It is shown that the agent's indifference prices 
are recursive if and only if the preferences are translation invariant. 
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Remark 5.1 Four definitions of g-risk measures are referred in this paper. They are 

<tE] UL 2 {F T -R n ); 

p 9 tiT [£] : = essinf {r, G L 2 (T t] R n )\ £j[^ + r ,\T t } > 0} , £ G L 2 (F T ;R n ); 
$ iT [£] := essinf {?? G I?(Tt\ R n )\ Sj[- (£ + t?) < 0} , £ G L 2 (F T \R n )- 

If g(t,y,z) = —g(t,—y,—z), then p 9 tT = ~p tT ; if g does not depend on y, then pf T = p 9 tT , 
PtT = PtT' if 9 does not depend on y and g(t,z) = —g(t, —z), then four definitions coincide with 
each other. Generally they are different definitions. We will continue our work to give insights 
among them and seek new applications in particular for the last two. 

6. CONCLUSION 

During a long time it seems that many problems related to multidimensional BSDEs are ex- 
pected to be complicated because of the lack of a clear structure of multidimensional BSDEs. This 
paper, by virtue of certain mathematical tools, establishes some properties of multidimensional 
g-risk measures. Particularly we show that a multidimensional g-risk measure is nonincreasingly 
convex if and only if the generator g satisfies a quasi-monotone increasingly convex condition. 

Indeed g-risk measures represent a large class of risk measures. In one dimensional case, any 
dynamic risk measure satisfying certain domination condition can be written as a conditional 
g-expectation (Coquet, Hu, Memin, and Peng (2002), Peng (2005)). As to the multidimensional 
case, similar result is presented in the appendix. Nevertheless, g-risk measures are interesting by 
themselves which are strongly linked to a functional g. Choosing carefully the coefficient g enables 
to generate dynamic g-risk measures that are locally compatible with the views and practice of 
different agents in the market. Especially a general dual representation of a multidimensional 
dynamic convex g-risk measure is obtained in which the penalty term is expressed explicitly via 
the polar function of g. It is shown that model uncertainty could lead to the convexity of risk 
measures. 

Our multidimensional approach provides an alternative way to measure the insolvency risk 
of a firm with interacted subsidiaries. As pointed out by Jarrow (2002), the cost of buying a 
put option written on the firm's net value is an intuitive measure of the insolvency risk of a 
firm. Optimal risk sharing between two groups with different risk tolerant coefficients is also 
investigated. We hope that this theory will be included in the toolbox of standard financial 
software in a near future. Different ways to obtain g-risk measures are compared at the end of 
the paper. 

Finally, there are also several interesting issues such as optimal risk transfer and risk allocation 
within a group with interacted subsidiaries, we do not consider due to the length of the paper. 
We will solve these issues in our future publications. 

Appendix A: BSDE REPRESENTATION FOR MULTIDIMENSIONAL INTERACTED 
RISK MEASURES AND NONLINEAR EXPECTATIONS 

This section will show that, in the setting of Brownian filtration, if a multidimensional interacted risk measure 
p Sj t satisfying certain domination condition, then it is in fact a g-risk measure, i.e. it can be represented by a 
multidimensional BSDE. For a multidimensional risk measure p s ,t = {pl,t, ■ ■ ■ , p",t) , we assume its components 
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are interacted with each other, or else it is just a notion putting n 1-dimensional risk measures together. Thus, for 
any k = 1, . . .,n, the feth component p k a is in fact a functional of the other dimensions, i.e., p k it depends on the 
whole path of (p') ; ^ fe on the time interval [s,i\. Precisely, we write p k t [-]( p i p k-i p k+i p n\ instead of p k tt [•] to 
put emphasis on the dependence of one dimension on the other dimensions. Sometimes for simplification we write 
an i? n_1 -valued vector (y t ) [ _ tk := (y 1 , ■ ■ ■ , y k ~ 1 ,y k+1 . . . , y n ). We make the following assumptions: 

(A6) Vfc = 1, . . . ,n, V (y l ), ,, G L^(0,r;i? n ~ 1 ), p k t [•]/ n is well defined on L 2 (T t ;R) and it satisfies (Al) 
~(A4). 

The above condition means that fixing other dimensions, the fcth dimension is still a risk measure satisfying 
(Al) ~(A4) except (A5). The condition of p 3,i -domination is 

(A7) For any k = 1, . . . ,n, 1 , ^ 6 L 2 (T t ; R), V (*J, . . . , a;?), {yl,..., yl 1 ) G L%(0, T; R n ), 




where p 9 a ^( k i ,\ is the solution of the following 1-dimensional BSDE: 

(A.l) Y a k = ~t + J^ M ( X)|»r| + |*?| + \Zr\ j dr ~ f s Z rdB r , s 6 [0, t]. 



Theorem A.l iet (Al) ~(A7) hold for the multidimensional dynamic risk measure p B ,t[']- Then there exists a 
function g : Q x [0,T] x R n x R nxd i — ► iT satisfying the "quasi-monotonicity condition" (3.4) and f/ie Lipschitz 
condition (H2) and g(, 0, 0) = 0, sttc/i f/iaf for any k = 1, . . . , n, 



We now present the proof in the language of nonlinear expectation. 

We define the nonlinear expectation operator £ s .t [£] := p s .t [—(,], which means that £ s . t satisfies (A'l) Mono- 
tonicity. Vt G [0,T], fi,,^ 1 ] > £ s>t [£ 2 ], if V? 1 > £ 2 € L 2 (T t ;R n ); (A/2) = £, V£ G L 2 (T t ;R n ); (A'3) Time- 

consistency. £., T [£] = £f, t [£ t B iT K]], VO < s < t < T, V£ G L 2 (T T ;R n ); (A'4) Regularity. U^.tKlx] = U£.,tK]> 
\M G J-^; (A'5) Normalization. £ s ,t[0] = 0; (A'6) Vfc = l,...,n, for any (n — l)-dimensional vector process 
(yl, . . .,y k -\y k+1 ...,y' s l )e L%(Q,T;R n ), £ k , t [•]( 1/ i,..., s }-i >1) J+i... > , ) ») is wdl defined on L 2 (T t ; R) and it satisfies 

(AT) ~(A'4); (A'7) f 9 " -domination. For any k = 1, . . . , n, V^ kl ,t 2 £ L 2 (F t ; R),"i (x\, . . . ,x n a ) , (yl, . . . ,y^) G 
2&(0,T; 7?™), 

(A.2) ^] -^ t [^] ^^[^-^] • 

where £ a JX [£ k ] , ,v is the solution of (A.l) with terminal datum £ fc G L 2 (Tt; R). 

The condition (A'5) means that £ Sjt [0]^ % = 0, which combining (A.2) yields that for any k = 1, . . . , n, 
G L 2 (F t ; R),\f(yl,..., y n s ) G L%'(0, T; R"), 



(a.3) ^[e*] 



1 (»o v * " • J (» ! U 

■WW*" ' ' ' 

(A.4) y s fe 



Let fC fe l 

/ ,\ , k — 1, . . . n, denotes the solution of the following 1-dimensional BSDE: 



t - [ fJ-\J2 \y'r\ + \ Y r\ + \Zr \ I * ~ / Z*dB T , S G [0,*]. 



Observing that £ s t **' [•]/ n = —£%1' [— \t n . Then by (A.2), one can directly check that: 
(A.5) -^] (a ,_ sf) <^ [^] fwiy 
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and 

(A.6) ^'' fe k fc ln <£s,tk k ] fl , ■ 

Let -Djr(0, T; R n ) denote all the processes in L 2 T (Q,T\ R n ) with right-continuous and with left-limit (RCLL) 
paths such that E [sup 0<t<T \>fit\ 2 ] < oo. 

Definition A.l Let E s .t be an n- dimensional nonlinear expectation. A process (Y s ) G 0^(0, T; R") is called 
£ -supermartmgale (£ -submartingale) if 

(A.7) £* t [Y t ] (Yl] < Y k , ("> " resp. ), Vs G [0, t] , h=l,. . . n. 

Obviously if £ s . t is £ s ^-dominated, then £ s . t is an £ s j -submartingale and an £ s t ^-supermartingale and £ s \ is 
an £ S!t -supermartingale. 

For a multidimensional nonlinear expectation operator £ and a process A 6 D%(0,T; R n ) and a random 
variable f £ L 2 (Tt\ R n ), we write £ t ,T [£; A\ := £ t ,r [£ + (At — At)] for brevity. To prove Theorem A.l, we need 
the following nonlinear Doob-Meyer decomposition. 

Proposition A.l Let (HI) ~(H3) hold for g and g k does not depend on (zj) ■ Let (Y t ) G D 2 T (0,T;R n ) be an 
£f T -supermartingale. Then there exists a unique increasing process (A s ) G D%(0,T; R n ) with Aq — such that 

(A.8) Y t = £l T [Y T ;A], 

i.e. for any k = 1, . . . ,n, 

rk -irk | / h ( _ v rrk \ i . / , k »k\ / rjh . 



k m j j-> 

z. dB s 



(A.9) Y t k = Y T R + J g K {s,Y s , Z K S ) ds + \A& - A$ J - J Z k dB 3 , 0<t<T, 

Proof. The main idea is to construct the following family of BSDEs parameterized by m=l,2, 
(A.10) V = Yt + f 9 k (s, Y s \..., fc y™ Y s k+1 Y?, k >z?) ds + m f ' (Y 3 k - k y?)ds - 



where k = l,...,n, to approximate the solution triple (Y,Z,A) of BSDE (A.9). Note that (A. 10) is in fact a 
1-dimensional BSDE. Then we can show that the sequence { k yT} _, is increasing in m and bounded by k y\ 
and Y k in L%(Q,T\ R) and k y™ converges up to the process (Y k ). The proof is an analogy of Lemma 3.4 and 
Theorem 3.3 in Peng (1999). We omit it here for brevity.. □ 

Corollary A.l For a given (C t ) G D%(0, T; R n ), if (Y t ) G D%(0, T; R") is an £ g t T [■;C]-supermartingale, i.e. 
(A.ll) k £l t W k -C k ] <Y k ,k=l,...n. 

Then there exists a unique increasing process (A s ) G Z)jr(0, T; _R") with Ao = such that 
(A. 12) Y t =£ 9 tT [Y T ;C + A]. 



Since (Ys' X '°) se ^ ( i is an £ft [•> C] -submartingale and an £ s t 9M [■; C*]-supermartingale, as a consequence of 
Corollary A.l and condition (A. 2), we have 

Lemma A.l Let (A'l) ~(A'7) hold. Given a process (C t ) G D%(0, T; R n ), for the process Y^ x - C := £ s>t [X;C], 
with s G [0,t], X G L (Ft; R n ), there exist a (gl' x ' c ) G L%(0,t;R n ) and a (Zl' x ' c ) G L%(0, t; R nxd ) such that, 
for any k = 1, ... ,71, 

(A.13) kyt,x,c = X k + ^ C k _ C k^ + j k g t,x,c dr _ j k z t,x,c dB ^ o < s < t, 
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and 

(A.14) |*fli ,x ' c |<ii(|Kf ,x,c | + |*^ ,jr,c |),Vfleto,t]. 

Moreover if (X, C ) changes, then 

j A -ir\ \k t.X.C k t,X' ,C ^ ( \ -\rt,X ,C -irt.X ,C i \k ryt.X.C k rjt.X' .Cl\\ w _ Xr\ j.1 

(A.15) \ 9b - 9s \<V\ys ~Y a | + | Zs ~ Z a ' J,Vs6[0,t]. 

We refer to Proposition 6.6 in Peng (2005) for similar arguments. 
Corollary A.2 Fort' e [t,T], X' g L 2 (T t >;R n ), (C' t ) G _D^(0, T; 7? n ), fe< y/.^'.c" := £ s t , [X';C]. T/ien 

Ifc t,X,C fe t'.X'.C'l ^ (\ v t,X,C v t'.X',C'\ , \kr,t.X,C k ryt' ,X' ,Cl\\ w . r„ ,, 

g e ~9s \<H\]Ys -Y a ' \ + \ Z„' ' - Z s M,Vs6[0,*J. 

Proof. It is just a consequence of the time consistency: Y s * ' x ' c = £ s ,t \y{ X ° ! C'j j s £ Pit] an d Lemma 
A.l. □ 

We now present a general nonlinear Doob-Meyer decomposition for £-supermartingale. 

Proposition A. 2 Let (A'l) ~(A'7) hold. Then an £-supermartingale (Y t ) G D 2 T (Q,T\R n ) has a unique decom- 
position 

(A.16) Y t = £ t , T [Y T ;A],t € [0,T], 

where (A s ) G Z?3r(0,T; R n ) is an increasing process with Aq = 0. 

The proof of the above proposition is quite similar to Proposition A.l and to Theorem 8.1 in Peng (2005). We 
need to construct the following sequence parameterized by m=l,2,. . . 



(A.17) k y? = £ 



t.T 



Y k +m l T (Y s k - k yT)ds 



(*2), 



The extension of £t,T to £ CTiT under stopping times is necessary. Then following the procedure in Theorem 8.1 in 
Peng (2005), the proof could be complete. 

Proof of Theorem A.l: For fixed (t, y, z) G [0, T] x R n x R nxd , consider the following n-dimensional SDE: 

n 

(A.18) d k Y^' z " = -^J2\ l Yl' y ' Zl \ + \z k \)ds + z k dB s ,se [t,T], 

l=i 

(A. 19) Y^' x " = y k , 

where k = 1, . . . , n, (yfc, Zk) € R x R d . By the classical theory of SDE, we have the following estimate: 

(A.20) E^yfYl"'"* -y k \ ] < {\y k \ 2 + \z k \ 2 + l) (s - t) ,Va G [t y T], 

fe=i fc=i 

where c is a universal constant changing from line to line in this section. 

By (A'7), (Ys' v ' Zk ) gS ^ t j is an £ s .t [-J-supermartingale. Then due to Proposition A.2, Y*' v '" k has the following 
decomposition: for any k = 1, . . . , n, 

(A.21) k Y t,y,z h =£sT \kyt,y,z k . A t, y ,z k ^ ^ g 

where A*' y ' Zfc G D%(0,T; R n ) is an increasing process with Ag V,Zk = 0. Applying Lemma A.l to (A.21), then 
there exists a pair ( k gl' y ' Zk , k zl' v ' Zk )£ L%(0,T;R) x L 2 T (0,T;R d ) such that 

(A.22) d k Y^ v ' Zk = -^'^ds-dA^ + k Z t a ' y ' Zk dB s , s G [t,T], 



34 



YUHONG XU 



and for each (y,z k ), {y',z' k )e R n x R d ,t<t' G [0,T], 

n 

(A.23) \ k gl' v ' Zk - k g i s' v ''*' k \ < fx^T \ l Y*' y > z > - <y/ V ' z '' | + ^>z* s ' y '*" - k zt'' v '' z ' k \),Vs e [t',T], 

1=1 

n 

(A.24) \ k gl' y " k \ < v{J2\' Y s' y ' Zk \ + |'^' y ' Zfc |),Vs G [t,T]. 

1=1 

For X G L 2 (T t ,;R n ), let F*'' X := £ s>t / [X] = £ Sjt , [X; 0]. Then by Lemma A.l and Corollary A.2 again, F*''* has 
the following decomposition, for any fc = 1, . . . , n, 

(A.25) d fc F*'' x = - k g^' x ds + k zt' x dB s , k Yt,' X = X k ,Vs G [0,*'] 

where the pair { k gi' x , k Z^ x )e L 2 T (0,t';R) x L 2 T (0, t'; R d ) and 

n 

(A.26) | fc ff^' Zfe ~ V S ' ,X | < ME l^'*' - lyt s' X \ + \ k Zl' y ' Zk - fe zl'' X |),Vs G [*,*'], 

i=i 

Comparing (A. 18) with (A. 22), we obtain h z\' v ' Zk = z k l[ t , T ] («)• Therefore (A.23), (A.24) and (A.26) could be 
rewritten as 

n 

(A.27) \ h gt' v ' Zk - h gl'' y ' A \ < m(E \' Y '' V ' Z ' ~ A \ + |** -4|),Vs € [t',T], 

Z = l 

n 

(A.28) | fc o^| <A*(£|'l?'H + |* fc |),V*e [t,T], 

1=1 

n 

(A.29) | fc ^' y ' Zfc - k fs' X \ < m(E \ l Ys' y ' Zl - l Yf' x \ + \z k - fc ^'' x |),Vs G [t,t'}. 

i=i 

For each natural number m > 1, we set if = i2~~ m T, i = 0, 1, 2, . . . 2 m , and define 

2 m -l 

(A.30) V" («, !/,**) == E k 9T' v '* k -l ltT , tT+l] (s), s€[0,T}. 

i=0 

Then k g m (s, y, z k ), k = 1, . . . , n, is J^-adapted and 

Lemma A. 2 For any fe = 1, . . . , n, /or eac/i /tred (j/, z^) G R" x i? d , t/ie sequence { k g m (■, y, «fc)} °_ 1 converges 
to some Ts-adapted process k g(-,y,z k ) G L^r(0,T;-R). 

Proof. For two natural numbers < m < m' , for each s G [0,T], there are some integers i < 2 m — 1 and 
j < 2 m ' - 1 such that s G n WP' ^7+i]- Then b y ( A - 27 ) an d (A. 20) we can deduce that 

sup 1? \ k g m (s, y, z k ) - k g m ' (s, y, z k )\ < 2ncf E (M* + \M 2 + l) h' m + 2' m ') T -> 0, 
36[0,T] I I fc=1 ^ / 

as m, m/ — > oo, which means that k g m (s, y, z k ") converges to some J^-adapted process k g (s, y, z k ) G L 2 T (Q, T; 7?). 

Lemma A. 3 For any k — l,...,n, the limit function k g(t,y,Zk) satisfies (i) \ k g (s,y, z k ) — k g {s,y , z' k )\ < 
»{\y-y'\ + \z k -z' k \), \fy,y' G R n , z k ,z' k G R d ; (ii) k g (s,y = 0, z k = 0) = 0; (Hi) \ k g (s,y, z k ) - k f X \ < 
v(\y-Y t s X \ + \ Zk - k zl' X \), Vs G [0,t]. 
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Proof. When s g [if, tj+i], by (A.30) and (A.27) we have 

n n 
fern/ \ fern/ / / \ ^ /I 'lil ' I \ i /\ ^ fcx/- £ 7 l 'f' 2 fc \ i /\ ^ 

3 (s,y,zk)- g [s,y ,z k )\ < fj,[\y-y \ + \z k - z k \) +^(> j Y s l - y k ) + n{) ^ 



k v t T' y '' z 'h l 



Then by (A. 20) we can show that the last two items of the right hand converge to zero as m — > oo. Thus (i) is 
obtained. Combining (A. 20) and (A.28), we could deduce (ii). Applying (A. 29) and (A.20) to j k g (s, y, z k ) - k g t a ' x \ , 
then (iii) is obtained by passing to the limit. See an earlier version of this paper in arxiv for complete proofs of 
the above two lemmas. 

Continuation of the proof of Theorem A.l: For each X G L 2 (J r t ; R n ), by (A. 25), F*' X := £ s , t [X] has the 
following decomposition: 

kT7*,X v k 



Y*,' X = X k + f tk g t ,: x dr+ [ k Z*; X dB r , s <E [0, t], fc = 1, 

J s J s 



By Lemma A. 3, the following n-dimensional BSDE, 



k Y*' x = X k + 



(r,Y,^ x dr - J k Z t r ' x dB r , s 6 [0, t], k = 1, . . . , n, 

t,X yt,X\ _ (l v t,X n\st,X 1 yt,X n yt,X\ 



has a unique pair of solution (F^ , ) = (T^ , . . . , n Yg , L Z\ 
to |F*' X - Y*' x \ , we get 

E^Yl' X - Y l s ' x ^ + E |z'.' X - Z*' x | 2 dr <2(/in + J E |F*' X - F r *' x | 2 dr + \e J ^ \zl' X - Z*' x \ 2 dr. 

It is a consequence of application of GronwalFs inequality that Y^ = Yg' x . Thus £ Stt [X] — Y f s ' X = Y^' x = 
t [X], where the function g satisfies (i) and (ii) in Lemma A. 3. Since p s . t [— ■ ] := £ s .t [•] is a quasi-monotone 

risk measure, inequality (3.4) holds also for g. □ 

Note that Theorem A.l could not be obtained directly from the one dimensional case. 
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